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ABSTRACT 

We present proper motions for the Large & Small Magellanic Clouds (LMC & SMC) based on 
three epochs of Hubble Space Telescope data, spanning a ^ 7 yr baseline, and centered on fields 
with background QSOs. The first two epochs, the subject of past analyses, were obtained with 
ACS/HRC, and have been reanalyzed here. The new third epoch with WFC3/UVIS increases the 
time baseline and provides better control of systematics. The three-epoch data yield proper motion 
random errors of only 1-2% per field. For the LMC this is sufficient to constrain the internal proper 
motion dynamics, as will be discussed in a separate paper. Here we focus on the implied center-of- 
mass proper motions: fJ,w,LMC = — 1.910 ± 0.020 mas yr~^ , ^at^lmc — 0.229 ± 0.047 mas yr~-^ , and 
Mw,SMC = —0.772 ± 0.063 mas yr~^ , ij,n,smc = — 1.117 ± 0.061 mas yr^^ . We combine the results 
with a revised understanding of the solar motion in the Milky Way to derive Galactocentric velocities: 
'ytot,LMC = 321 ± 24 km s~^ and ?;tot,SMC = 217 ± 26 km . Our proper motion uncertainties are 
now dominated by limitations in our understanding of the internal kinematics and geometry of the 
Clouds, and our velocity uncertainties are dominated by distance errors. Orbit calculations for the 
Clouds around the Milky Way allow a range of orbital periods, depending on the uncertain masses 
of the Milky Way and LMC. Periods < 4 Gyr are ruled out, which poses a challenge for traditional 
Magellanic Stream models. First-infall orbits are preferred (as supported by other arguments as well) 
if one imposes the requirement that the LMC and SMC must have been a bound pair for at least 
several Gyr. 

Subject headings: galaxies: interactions — galaxies: kinematics and dynamics — galaxies: evolution 
— Galaxy: structure — Magellanic Clouds 
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1. INTRODUCTION 

High-precision proper motion (PM) measurements of 
the Large and Small Magellanic Clouds made by our 
group with two epochs of Hubble Space Telescope's 
(HST) ACS High Reso lution Camera (HRC) data 
(|Kallivavalil et all l2006aira . hereafter K2 & Kl), and 
confirmed bv iPiatek etalf ()2008l . hereafter P08), have 
revolutionized the field of Magellanic Clouds research. 
The implied tangential velocities were high enough, ap- 
proximately 100 km s~^ higher than in previous theo- 
retical models, and the observational errors were small 
enough 0.07 mas yr~^ for LMC), that an orbital so- 
lution in which the Clouds are either only now on their 
first infall about the Milky Way (MW) ()Besla et al.ll2007l 
hereafter B07), or arc on an eccentric, long period (> 6 
Gyr) orbit ( Shattow &: Loeb 2009, B07) are the favored 
solutions. In addition, the observed relative velocity be- 
tween the Clouds was of order the escape speed of the 
SMC from the LMC (105 ± 42 km s'^; see K2) leaving 
open the possibility that the Clouds may not be bound 
to each other, although bound orbits were still allowed 
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within the relatively large error bars which came from the 
less-precise SMC PM determination (~ 0.18 mas yr~^ ). 
These results received much attention from the commu- 
nity in part because they require a new formation mecha- 
nism for the Magellanic Stream, a young coherent stream 
of H I gas tha t trails the Clouds ^ 150° across the 
skv (IWannier fc Wrixon 1972; Mathcwson e tHI 119741: 
Putman et al. r T2003t iBriins et al.. .2005^ .Nidever et al l 
2010| ). Most models, be they tidal or ram-pressure in 
nature, require multiple pericentric passages in order 
to b e viable s tripping mecha nisms ([Rtizicka et al. 2003 
Bekki 2008; Connors e t al. 2006: Mastropietro et 3] 
2005: Gardiner fc Noguchl .1996: Moore fc Davis. ,1991 
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We have recently put forth an alternative formation 
mechanism for the Stream in which the material is re- 
moved by LMC tides acting on the SMC before t he sys- 
tem falls into the MW for the first time (Besl a et al.l 
120101 ). A firm prediction of this model is that the Clouds 
have been bound to each other and further that the SMC 
is on a highly eccentric, prograde orbit about the LMC 
(e — 0.7; apocenter ~ 100 kpc). This orbital configura- 
tion prevents the Clouds from merging and also leaves the 
dispersion-supported material within the SMC's disk ra- 
dius relatively unaffected, while resonances aid in remov- 
ing rotationally supported material from the outskirts of 
the SMC's disk. Thus, despite the large parameter space 
involved in modeling a feature such as the Stream, there 
are two pieces of information that can dramatically re- 
duce the uncertainty in the models: 1) knowing whether 
the Clouds are on a short or long period orbit about 
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the MW, and 2) the orbital eccentricity of the SMC's 
orbit about the LMC. We rece ntly expanded on this 
work, building on earlier ideas o f lYoshizawa &: Noguchil 
12003) and Be kki fc Chibal &K)Fh . arguing that the inter- 
nal kinematics and structure of the LMC strongly favor 
a scenario in w hich the MCs have recently experienced a 
direct collision (jBesla et al.ll201^ . We also looked more 
generally at the implications of such dwarf-dwarf interac- 
tions for Magellanic Irregulars, a class of dwarf galaxies 
for which the LMC is a prototype. We found that pro- 
grade dwarf-dwarf tidal interactions can efficiently re- 
move baryons from the lower mass companion and that 
structures such as off-centered bars and one-armed spi- 
rals may be hallmarks of ongoing or recent interactions 
with a low mass companion (mass ratio 1:10). 

Apart from our own work on the formation of the 
Stream, these PM measurements have spurred many 
other new theoretical investigations into the Clouds' ori- 
gin, wh ether from M31 or the far reaches of the Local 
Group (iShattow fc Loeb| [20Ql iKallivavalil etall l2009al: 
lYang fc: Hammed 120101 ). Attempts have been made to 
measure the Milky Way mass from the kinematic prop- 
erties of the LMC (iKallivaYahl et al. 2009a; Busha et al.l 
l2011at iBoylan-Kolchin et al J 2011t ) . Recent studies have 
attempted to place the orbits of the Magellanic Clouds 
in a cosmological context, looking at the expected infall 
times in a statistical sense, and favoring a scenario in 
which the Clouds are likely on their first infa,ll into the 
MW (Bovlan-Kolchin et al. 2011:'Bu sha et al.l[20TTR see 
also Bek ki (2011.) ). Rocha et a l. (2 0lll ) analvzed subha- 
los in the Via Lactea II cosmological simulation and find 
that present day orbital energies are tightly correlated 
with the time at which subhalos crossed into the host's 
virial radius. Their analysis indicates that the LMC en- 
tered th e MW virial r adius ~ 4 Gyr ago, in agreement 
with t he iBusha et al.l (|2011al ) and iBovlan-Kolchin et al.l 
(|20Tll) results. 

Overall, the IBusha et al] (|2011aD and 
IBovlan-Kolchin et al.l (|201lD works find different 
conclusions for the mass of the MW. This is because 
Busha et al. also folded in constraints on the separation 
between the Clouds and the MW host which yields 
posterior distributions that favor lower mass MW halos. 
The Boylan-Kochin et al. study favors MW models with 
a high mass because the Kl velocities imply the LMC 
is moving close to the escape speed of a IO^^Mq halo. 
Since subhalos are rarely ever accreted on unbound 
orbits, this means that the Bayesian analysis will be 
weighted towards a higher mass. Consequently, they 
favor a MW mass in excess of 2 x 10 ^^M©. Both 
studies represent a novel approach and new avenue 
for near- field cosmology, i.e., using accurate orbital 
kinematics/histories of MW satellites to constrain the 
properties of the MW. Because the speed is so close 
to the escape speed it is crucial that the velocity 
measurement be accurate for the LMC. A lower speed 
for the LMC would favor MW mass estimates between 
1-2 xIQI^Mq. 

In order to better address these questions about the 
Clouds' velocities and likely orbits, we obtained an addi- 
tional epoch of HST data with the Wide Field Camera 
3 (WFC3). With measurement errors similar to those 
in epochs 1 & 2, and an increased time baseline from 2 
to 7 yrs, this allows for decreased PM errors by a factor 



of a few. In turn, this yields improved knowledge of the 
past orbit. Moreover, a third epoch provides a valuable 
check on systematic errors: there is always a straight line 
between two points, but if a third point, particularly ob- 
tained with a completely different detector, fails to line 
up then that is a clear sign of errors in the analysis. 

In this paper we present the results from the analysis 
of the third epoch of HST data of the Magellanic Clouds, 
and the implications for their space motion and past or- 
bit. In f|2]we discuss the details of the observations, the 
analysis of the WFC3 data, an improved reanalysis of 
the ACS data, and the resulting two-epoch and three- 
epoch PMs. In P we derive the PMs of the LMC and 
SMC Center of Mass (COM) from the measurements of 
the individual fields. In 21 we compare the new results 
to previously reported measurements and estimates. In 
f|5] we discuss the corresponding space motions of the 
Clouds. In Sj6]we discuss the implications for the orbits 
of the Clouds, and the new insights and improvements 
thus obtained. The main conclusions are summarized 
and discussed in fjT] Appendices discuss some additional 
HST data that we obtained, but which are not included 
in the present study. 

In a companion Paper II we use the new LMC data to 
study its PM rotation field. This provides new insights 
into the distance, center, orientation, and rotation of the 
LMC disk. The COM motions derived in the present 
paper use the field dependent corrections derived in Pa- 
per II. 

2. DATA & ANALYSIS 
2.1. Description of Observations 

iGeha et al.l (|2003( ) identified a total of 54 QSOs behind 
the Clouds from their optical variability in the MACHO 
database. These QSOs provided the inertial reference 
frame against which the PM measurements were made 
over a 2 yr baseline in Kl & K2, based on observations 
obtained between August 2002 and June 2005. For most 
efficient use of HST resources we decided to observe in 
SNAPSHOT mode for both epochs. Our two SNAP pro- 
grams of imaging with the HRC achieved an overall com- 
pletion rate of 48% yielding a final sample of 21 QSOs 
behind the LMC and 5 behind the SMC. 

We obtained an additional epoch of data during the 
period July 2007 to November 2008. These observations 
executed with the HST WFPC2 camera, due to the fail- 
ure of ACS at that time. We briefly describe these data 
in Appendix A, but due to their limited quality, we don't 
use them in the analysis presented here. 

We obtained another epoch of data with the HST 
WFC3 camera during the period October 2009 to July 
2010. These observations also executed in SNAPSHOT 
mode, and we refer to them as the "third epoch" . The 
final yield from this program was 11 observed QSO fields 
for the LMC and 4 observed QSO fields for the SMC. 
However, two QSO fields, one in each galaxy, proved 
unsuitable for PM determinations. As discussed in Ap- 
pendix B, in one field the QSO showed a bright extended 
host galaxy, which complicates astrometry, and in the 
other field the nature of the QSO was not confirmed. 
These two fields are omitted from the subsequent discus- 
sion and analysis. Therefore, the final number of QSOs 
used for the three-epoch analysis is 10 for the LMC and 
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3 for the SMC. 

The HRC data, which makes up the first two epochs 
of this study, is discussed in detail in Kl. Briefly, for 
our main astrometry goals we chose the F606W filter 
which is a broad V filter with high throughput. Eight 
dithered exposures were taken in each epoch to minimize 
pixelization-related systematic errors, and the dither pat- 
tern was kept the same over the two epochs. Exposure 
times were chosen so as to achieve a S/N of at least 100 
for the QSOs based on the known MACHO magnitudes 
at the time. 

For the third epoch with WFC3/UVIS (PID 11730), 
we again used the F606W filter and used the 4-point 
DITHER-BOX pattern which is made for optimal sam- 
phng of the PSF. We aimed at S/N ~ 200 for the 
QSOs in order to match the astrometric quality of the 
ACS data. Even though the WFC3/UVIS pixels are big- 
ger than ACS/HRC pixels, with adequate dithering this 
should not degrade the astrometry because the astromet- 
ric error is proportional to the [FWHM/ SNR) of the 
target. With QSO brightnesses in the epoch 1 sample 
ranging from 16.5 <V< 22.0, and an average brightness 
oiV — 19.5, this yielded total science exposure times for 
epoch 3 ranging from 2.6 to 17.7 minutes, with an aver- 
age of 12 minutes. Tables [T] & [2] describe the combined 
three-epoch dataset, including RA/DECs for the QSOs, 
whether they were observed with WFC3 or not, the to- 
tal WFC3 exposure times, and the corresponding time 
baseline for the PM measurement. 

Figure [T] shows the locations of the QSOs behind the 
LMC and SMC. QSOs for which we obtained two epochs 
of data are marked with yellow diamonds. Their distri- 
bution behind the LMC is sparse but reasonably uniform. 
QSOs observed with the WFC3 are marked with red 
squares. Note that there is 1 "new" LMC field for which 
there is only first epoch ACS and third epoch WFC3 
data. 

2.2. Analysis of WFC3/UVIS Data 

As in Kl & K2, we used the bias-subtracted, dark- 
subtracted, flat-fielded images (_flt.fits) provided by the 
STScI da ta reduction pipel i ne for our basic reduction 
purposes. lAnderson fc Kind ()20Q4[ ) found that these im- 
ages are well-behaved astrometrically. They are not, 
however, corrected for geometric distortion. We did not 
use the geometrically corrected products created by the 
MultiDrizzle software (_drz.fits) since they involve a re- 
sampling that might degrade the astrometry. Instead 
of geometrically correcting the images, we geometrically 
corrected the positions that we measured on the images. 
The positions were measured by means of empirical PSFs 
constructed from other data sets (similar to the approach 
discussed by Anderson & King 2006 in ACS/ISR 06-01). 
The geometr i c corr ections were performed according to 
iBellini et al.l (|201lD . In their study, the geometric dis- 
tortion is modeled with two components: a third-order 
polynomial and an empirically-derived look-up table for 
the residuals, one per chip and filter combi nation. The 
method is akin to that used m,A nderson fc King! (|2004l ) 
and is shown in Be llini et al.l (|201ll l to be stable over time 
to better than 0.008 pixels (RMS) in each coordinate. 

Since we are performing relative astrometry we have 
freedom in constructing the masterframe (per QSO field) 
into which all relevant star positions are to be trans- 



formed. While PSF-fitting was performed on the fit im- 
ages, we did use stars in the drz image as a starting point 
to defining our masterframe, because their WCS header 
information is the most accurate. Simple aperture pho- 
tometry was performed on the drz image to obtain initial 
positions and magnitudes for real sources in the region 
of overlap of all 4 individual exposures with orientation 
North up and East to the left. We then did a sim- 
ple transformation of these initial positions to an (a;, y)- 
frame of identical orientation, but with a chosen center 
and 25 mas pixels (for simplicity). A six-parameter linear 
transformation for the 4 individual fit WFC3 exposures 
into this frame was then calculated. Based on previous 
experience, we only use stars brighter than an instru- 
mental magnitude, A/instr = -8 {S/N ^ 40). The 
instrumental magnitude is defined as — 2.5 log(coMnts) 
in the F606W fit image (without application of a zero- 
point). We then calculated a new masterframe position 
for each source by averaging the linearly transformed 
WFC3 positions. At this stage we iterated on the lin- 
ear transformations, using only stars with measurement 
errors smaller than ~ 0.1 pixels, that are common to all 
4 WFC3 exposures. The positional errors are calculated 
as the RMS scatter between multiple measurements of 
the same source. 

In general, we can routinely center WFC3 stars to 0.02 
pixels. Because our exposures are deeper and the FOV is 
much larger than for ACS/HRC, we identify on average 
~ 200 — 800 sources in the WFC3 images, depending on 
the field. By contrast there are far fewer sources in the 
HRC fields. 

2.3. Re-analysis of ACS/HRC Data 

In our original Kl analysis we argued, based on de- 
tailed calculations and consistency checks, that the de- 
grading charge transfer efficiency (CTE) of the HRC 
was not expected to systematically affect our analysis. 
The main reason was that the SNAP nature of the pro- 
gram yielded random roll angles of the telescope for each 
QSO field. Given that we had TV = 21 QSO fields with 
more or less random detector orientations on the sky, 
any CTE induced astrometric shift along columns would 
average down to zero oc iV-i/^ There is no explicit 
model for the underlying HRC CCD detector physics 
and char ge trap properties, as th ere now is for other 
cameras (| Anderson fc BedinI IMTol ) . The residual CTE 
effects, given the random roll angles of our exposures, 
would be expected to puff up the final RMS of the mea- 
surements. We included this effect by using the RMS 
scatter between fields as our final error estimates, in- 
stead of simply propagating the per-field weighted errors 
(which would have suggested a factor of 2 higher accu- 
racy). 

In Kl, our final PM value was based on 13 fields that 
we deemed as high quality for various reasons, and be- 
cause of this reduced number of fields, and the small 
number of total SMC fields, CTE issues remained a le- 
gitimate concern. P08 subsequently performed an inde- 
pendent re-analysis of our HRC data and used a sim- 
ple correction for CTE that was a function of Y co- 
ordinate of the object, the time since installation, and 
oc S/N~^-'^'^. This latter exponent came from the ACS 
Handbook. They found, as expected, that correcting for 
CTE did not have a large systematic effect on the anal- 
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Figure 1. i?-band image of the LMC (8° X 8°), and the SMC (3° X 5°). The MACHO photometric coverage is indicated. Red squares 
indicate reference QSOs with WFC3 data and yellow diamonds indicate those with two epochs of ACS data. 

ysis (they obtained a COM PM for the LMC within la 

from that of Kl). However, they did obtain better field- 

to-ficld agreement. This yielded smaller error bars on the : ' 

COM motion, and the possibility to derive a rudimentary : 

PM rotation curve. * : 

For the present study we reanalyzed the HRC data us- : 

ing our own new prescription for CTE degradation as : 

well as other improvements. The main differences in this a - 

reanalysis versus the analysis in Kl & K2 are: (1) we § 

use new codes for the implementation of 6-parameter lin- > : 

ear transformations between the two epochs that employ ^ '-_ 

error-weighting for the stars. In the previous analysis '- 

we only used rejection. (2) We employ a prescription to - 

correct for CTE, the form for which is as follows: : 



2;new = J/RAW - (a X 



y 



t 



1024.0 365.25 

MAX[(MiNSTR - MlNSTRJim), 0] . 



norm 



(1) 



where y is the measured position of a star on the de- 
tector, MiNSTR is the instrumental magnitude, typically 
ranging from Mjnstr = —6 (faintest) to —12 (bright- 
est) in our data, and t is the time since installation of 
ACS (Modified Julian Date = 52333). CTE affects faint 
stars more than bright stars, and stars brighter than 
-^iNSTR,iim are assumed to be unaffected. Norm is sim- 
ply a normalizing factor discussed below. The slope, a, 
was fit in the following fashion: we sequentially corrected 
averaged detector y positions by a range of slopes. Af- 
ter application of the correction of a given slope, error- 
weighted transformations were calculated for stars from 
epoch 2 into the epoch 1 frame. The input error per 
epoch is the RMS / ^/N obtained from centering the star 
over N exposures in that epoch. We then calculated what 



Figure 2. The total ^ jDOF for all the QSO fields as a func- 
tion of CTE slope Q (see Equation [TJ , for AfiNSTH,,iim = —10 and 
norm, = 4. The minimum indicates the value of a that provides 
the best correction. 



slope gave the minimum total / DOF for these trans- 
formations summed over all 21 QSO fields, where DOF 
stands for degrees of freedom, and is defined as follows: 
A^DOF = A^data " A^param, whcre A^data IS twicc the num- 
ber of fields and A^param is the number of free parameters 
optimized in the fit. The QSO is not included in these 
transformations . 

Figure [H shows the total I DOF for all the QSO 
fields as a function of a. We tried several MiNSTR.iim 
hmits, MiNSTR,iim = (—10,-11,-12), where —12 corre- 
sponds roughly to the brightest stars in the field. We 
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Figure 3. Residual star PMs (pix) after CTE correction via Equa- 
tion |[T} as a function of MmjgxR (black symbols) and plotted in 
the ACS el frame, i.e. x and y are in the native detector direc- 
tions. The pixel scale is 28 mas/pix. Stars from all fields were 
binned together. The average PM should be zero by construction. 
Remaining trends are indicative of systematic errors. Points in red 
show the results without application of CTE corrections. 

normahzed the magnitude correction such that it is the 
greater of zero, when a star is brighter than or equal to 
— 12, or a fractional correction up to a size of approx- 
imately 1, for AfiNSTR = —6 to —12. Therefore, for 
A'/iNSTR,iim = (—10, —11, —12), the Corresponding values 
for norm are (4, 5, 6). This range of MiNSTR,iim did not 
affect the final PM values by more than la. In addi- 
tion to fitting for CTE-slope, we tried a number of dif- 
ferent faint end magnitude cut-offs for the stars used in 
the transformations, settling on only using stars brighter 
than MiNSTR = -8. We also employed rejection in the 
form of (I residual/error I < 3), i.e. a cut based on the size 
of the residual of the transformation of a star divided by 
its measurement error. 

As a consistency check, we calculated the residual in 
pixels for each Magellanic Cloud star between epoch 2 
and epoch 1. This PM should be zero on average, since 
differential streaming motions within each field are too 
small to be resolved by our measurements. We averaged 
the y-residuals over all stars in all QSO fields, and binned 
the results by Mjnstr- The QSO was again excluded. 
Figure [3] shows the results. Points in red correspond to 
a = (i.e., no CTE-correction applied). A clear trend 
with magnitude in the y-residuals reveals CTE-induced 
astrometric shifts. By contrast the black points are for 
our best-fit CTE-correction parameters in eq. ([1]), which 
are a = 0.048 for MiNSTR,iim = — 10 and norm — 4. 
There is no evidence for any remaining y-CTE trend 
larger than 0.003 pixels (i.e., smaller than our random 
errors). There is no indication that serial-transfer CTE 
in the x-direction might have significantly affected the 
analysis, so we do not apply any corrections for it. 

2.4. Two-epoch Proper Motions 

From hereon we define the proper motion in the west 
{fj.w) and north {(In) directions in terms of the change 
in right ascension (a) and declination (6) on the sky: 



= — {da / dt) cos{6) , fiN = d6/dt. With the final ACS 
to ACS transformation terms in hand, we determined the 
PM for each source, including the QSO. A linear func- 
tion with time was fit to the available positions. The 
best-fit slope gives the PM. The average PM of the stars 
in each field is zero by construction. However, the QSO 
appears to move due to its reflex motion with respect to 
the LMC/SMC stars in the foreground. Hence, the aver- 
age absolute PM of the stars in the field is simply minus 
the PM of the QSO. The final error in this PM has two 
components. One component is the error in the average 
PM of aU the stars: cr(PM)0. This error quantifies how 
accurately we were able to align the star-fields between 
the epochs. This was added in quadrature to the for- 
mal error of the fit to the QSO's motion, ^PMqsq. So 
the error of a given field in either direction is 6fj,w or 

S^^N = ^JsPMqso^ + <j1pm)- 

Columns 7-10 of Tables [T] & H list the PMs thus in- 
ferred for all fields with two epochs of ACS data, and 
the associated errors, {Sfiw, ^Mat)- These results are very 
consistent with those of P08 even though we have used 
different methodologies for PSF-fitting, linear transfor- 
mations and CTE-correction. The average residual be- 
tween our reanalysis and that of P08 is 0.08 mas yr~^ 
lower in the West direction and 0.06 mas yr~^ lower in 
the North direction, while the summed over all 21 
fields is 12 in the West direction and 8 in the North 
direction. This is acceptable compared to iVnoF- The 
1(7 difference between our final values and those of P08 
arises from the inclusion of error- weighting in the trans- 
formations, and not the specifics of the CTE correction. 
If we exclude error-weighting, we recover values closer to 
P08. Since there are some sparse fields in which only a 
small number of stars are used in the transformations, it 
is especially important to include error-weighting in the 
analysis. 

2.5. Three- epoch Proper Motions 

The inclusion of the WFC3 data gives us 20 images per 
QSO field for 15 fields: 8 ACS first epoch (el), 8 ACS 
second epoch (e2), and 4 WFC3 third epoch (e3). One 
LMC field has only el and e3 data, and hence 12 images. 
To calculate PMs we first correct the ACS positions for 
CTE as in i i2.3l We then perform preliminary transfor- 
mations for stars (brighter than AfiNSTR — —8) in each 
of the 20 images into the masterframe defined in H2.2I All 
20 resulting {x, y) centers of each star are averaged and 
the RMS is calculated to make up the final masterframe 
and error (the QSO positions are not averaged) . We then 
treat this as the initial target frame for transformations. 
As the input frame, we use the individual positions of the 
stars in each of the 20 images, with input error the RMS 
per epoch. Again, as in ^2.31 error- weighting, rejection, 
and iteration are used in the determination of the final 
transformation terms into the masterframe. 

At this stage we lose the possible advantages of hav- 
ing such a large FOV and high number of sources in 
e3 because we are limited by the relatively sparse (and 
small) HRC fields. When we perform final transforma- 

^ Although the average is formally zero by construction, the 
precision of this zeropoint is limited by the measurement errors of 
the stars at the two epochs. 
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tions to put el, e2 and e3 into the masterframe, we only 
use the sources that are common to the HRC. The num- 
ber of sources used in the transformations, after all cuts, 
is listed in column 6 of Tables [1] & [2] The variation in 
this number reflects the real variation in stellar density 
at different locations in the Clouds and is not an artifact 
of the analysis. 

Figure |4] shows the x vs. y (equivalent to W, N) po- 
sitions in pixels of the QSO, after transformation into 
the masterframe, for 4 randomly chosen QSO fields. The 
fields are labeled according to Tables [T] & [51 Linear mo- 
tion is clearly visible, as is the consistency across the 
three epochs. Each triangle represents 1 of the 20 posi- 
tions of the QSO relative to the starfield. The scatter re- 
flects the real error in centering the QSO within a given 
epoch. In most cases the quality (RMS) of the WFC3 
data is comparable to that of HRC as we expected based 
on the fact that we aimed for higher S/N in e3 to account 
for the larger pixels of the WFC3. This plot represents 
the reflex motion of the L/SMC. The actual PM and its 
random error were determined as described in '> \2A\ by 
fitting linear motion with time to all sources. The high 
quality of the results is visible by eye in Figure 01 For 
two of the LMC fields in Figure SI we also show the x and 
y position in pixels of the QSO as a function of time in 
Figure [3 The corresponding linear fits are also shown. 

Columns 11-14 of Tables [T] & [2] list the PMs thus in- 
ferred for all fields with both ACS and WFC3 data, and 
the associated errors, {5^w i^^-n)- 

3. PROPER MOTION RESULTS 

3.1. Comparison of Two-Epoch and Three-epoch PMs 

In Figures m and [3( top) we show the measured QSO 
motions (colored symbols) in the (^vf, MJv)-plane com- 
pared to the star motions (open circles; centered around 
zero by construction) for the LMC & SMC respectively. 
The stellar motions clearly separate from the QSO mo- 
tions. The green filled circles show the results from linear 
fits to all three epochs of data, and the pink squares show 
results from ^2.31 for el and e2 only. The fact that the 
averages are consistent with each other is visible by eye. 
The addition of a third epoch with a different HST in- 
strument not only improves the accuracy of the measured 
PMs, but also demonstrates that there were no funda- 
mental systematic problems with the earlier two-epoch 
analysis of the observations in Kl, K2, and P08. 

3.2. LMC COM Motion 

As discussed in Kl, given the large size the LMC sub- 
tends on the sky and the fact that we know it pos- 
sesses internal rotation, the motion observed for an in- 
dividual QSO field must be written as: PM(ficld) = 
PM(COM)-}-PMi.es (field). Here PM(COM) is the proper 
motion of the LMC COM and PM^es (field) is a field- 
dependent residual. The latter contains contributions 
from different effects: (a) variations as a function of po- 
sition in the components of the three-dimensional veloc- 
ity vector of the COM that are seen in the plane of the 
sky ("viewing perspective"), and ( b) the internal rota- 
tion of the LMC (see formulae in Ivan der Marel et al.l 
120021 hereafter vdM02). Both corrections make a simi- 
lar contribution to the observed PMs. Given the loca- 
tion of our fields, the viewing perspective contributes 



< 0.4 mas yr~^ and the internal LMC rotation con- 
tributes < 0.3 mas yr^"'^ . This is less than ^ 20% of 
the overall PM, but that is significant given the accu- 
racy that we are trying to achieve. 

The field-dependent residuals depend on the COM PM 
and vice-versa and hence we used an iterative proce- 
dure to evaluate them in Kl. For the internal rotation 
model of the LMC we used the model of vdM02, based 
on the line-of-sight (LOS) velocity data for 1041 carbon 
stars. This was because our PMs themselves were not 
accurate enough to significantly constrain the LMC ro- 
tation model. P08 took a similar approach by leaving 
the LMC center and disk plane orientation fixed at the 
values found by vdM02. However, they did vary and fit 
the amplitude of the rotation curve. 

With the addition of the WFC3 epoch, the per- 
ficld, per coordinate PM errors are now typically 
0.03 mas yr~^ (see Table [U columns 13 & 14), which 
is only ^ 7 km s^^ . We show in Paper II that this pro- 
vides enough information to independently constrain all 
parameters of the PM rotation field, without reference to 
prior knowledge from LOS velocity studies. Table [31 lists 
the results, and compares them to the values from vdM02 
which were used by Kl and P08. We also co mpare them 
to the values from a more recent study by lOlsen et al.l 
(|20lH ). who analyzed the LOS velocities of ~ 6000 mas- 
sive red supergiants, oxygen-rich and carbon-rich AGB 
stars. A detailed discussion of the commonalities and 
differences between the results is presented in Paper II. 
We summarize here only the salient features relevant to 
the determination of the LMC COM motion. 

The inferred orientation of the LMC disk, as mea- 
sured by the inclination i and the position angle 9 of 
the line of nodes, is consistent with the range of values 
previously rep o rted u sing other methods (see reviews of 
Ivan der Marell (|2006t ) and Ivan der Marel et al.l (|2009[) V 
The rotation curve amplitude inferred from the PM data 
falls between those determin ed from the LOS velo cities 
of (old) carbon stars (vdM02: Olson fc Massev ll2007D and 
(young) red supergiants (Olson ot al. 2011, see Table 3). 
This is as expected, since the (bright) stars in our HST 
fields contain a mix of older and younger stars (see the 
color- magnitude diagram in Figure 6 of Kl). This con- 
firms the accuracy of the new PM data, both in a random 
and a systematic sense. The rotation curve amplitude in- 
ferred from our PM data is not, however, consistent with 
the result of P08. They were able to infer a rotation 
curve from the two-epoch data, but their rotation ve- 
locity was surprisingly high: 120 ± 15 km . This is 
approximately 30-40 km s~^ higher than the value de- 
rived from the LOS velocities of H I and red supergiants 
()Kim et al . 1998; Olscn ct al. 2011). It would be particu- 
larly hard to understand how the stars in the LMC could 
be rotating significantly faster than the H I gas. 

The most important difference between our new re- 
sults from Paper II and the LOS study of vdM02 (used 
by Kl and P08) is in the position of the LMC center. 
In Paper II we find strong evidence that the center of 
the stellar PM rotation field is consistent with the po- 
sition of the H I dynamical center. This makes theo- 
retical sense, since the stars and gas are orbiting in the 
same gravitational potential. However, this result differs 
by 1.12° from the position advocated by vdM02, which 
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Figure 4. The x vs. y (equivalent to W, N) positions in pixels of the QSO, after transformation into the masterframe, for 4 randomly 
chosen QSO fields for all 3 epochs of data. The fields are labeled according to Table [l] Linear motion (reflex motion of the L/SMC) is 
clearly visible, as is the high S/N consistency across the three epochs. Each symbol represents 1 of the 20 positions of the QSO relative to 
the starfield, and are distinguished here using triangles for ACS data and plus signs for WFC3 data. The uncertainty in centering the QSO 
within a given epoch is visible as scatter between the points. The first two ACS epochs are the ones that are most closely spaced in time 
(and QSO position). The QSO's reflex motion is predominantly towards the West, because the Magellanic Clouds move predominantly 
towards the East. Each panel is 1 X 1 pixel. 1 pixel = 25 mas. 



agrees with the brightest part of the LMC bar. The use of 
a new rotation center affects the inferred PM of the LMC 
COM, {^iw,LMCi l^N ,LMc), as discussed in Section HTTl 

In Table [T] we list for each field the estimate for the 
LMC COM that is obtained after subtraction of the con- 
tributions from viewing perspective and LMC rotation 
indicated by the best fit model. In Figure ^[bottom) 
we show the estimates PMost(COM) resulting from this 
procedure. The scatter is considerably reduced with re- 
spect to the top panel. The weighted average of the 
estimates is hwxmc — —1.910 ± 0.008 mas yr~^ and 
fJ-N.LMC = 0.229 ± 0.008 mas yr~^ . The error bars on 
these values reflect only the propagation of the random 
errors in the data. We show in Paper II that there are 
also contributions from two other sources. First, there 
is excess scatter between measurements from different 
fields that is not accounted for by random errors, disk ro- 



tation, and viewing perspective. The best-fitting model 
has a reduced (XminZ-^DOp)"^^^ — 1.8; here Xmin is the 
value of the best-fit model. Second, there are uncer- 
tainties in the geometry and rotation of the best-fitting 
LMC disk model that need to be propagated. The lat- 
ter provide the dominant uncertainty in the final result. 
When these sources of uncertainty are taken into account 
as described in Paper II, the final estimate of the LMC 
COM PM becomes ^lw,LMC = -1.910± 0.020 mas yr"! 
and fiNXMC = 0.229 ± 0.047 mas yr^^ . 

3.3. SMC COM Motion 

For the SMC we have PM measurements for only 5 
QSO fields, of which only 3 fields have three-epoch data. 
This sparse coverage cannot provide much insight into 
the SMC geometry and rotation field. We therefore fit 
a relatively simple model to the PM data in order to 
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Figure 5. Analogous to Figure|4]but now showing x and y positions of the QSO in pixels as a function of time for two of the LMC fields 
also in Figure|4] The fitted linear motion is depicted by the solid lines. 1 pix = 25 mas. 



determine the SMC COM PM {fiw^ 

SMC, fJ-N,SMC)- 

To calculate a PM model prediction for each observed 
field, we keep the SMC center fixed at t he H I kinemati- 
cal ce nter (a, S) = (16.25°, -72.42°) (St animirovic et al.l 
120041) ■ the distan ce modulus fixed at m - M = 18.99 
(jCioni et al.ll2000l) , and the radial velocity fixed at Vgys = 
145.6 km (|Harris fc Zaritskvl[2006l ). We account for 
the influence of viewing perspective as in vdM02. We 
also allow for the possibility of a single overall rotation 
of velocity V^ot in the plane of the sky (i.e., as though we 
were viewing a face-on disk). We treat the SMC COM 
PM as a free parameter that is optimized by minimizing 
the of the model fit to the data. 

To determine the uncertainties on the best-fit model 
parameters, we create pseudo-data in a Monte-Carlo 
sense from the best-fitting model, with properties sim- 
ilar to the real data. For this we use the observational 
error bars, but multiplied by a factor {Xmin/'^DOFy^^ 
to account for the observed scatter between the mea- 
surements. Many different pseudo-data sets are created 
that are analyzed similarly to the real data set. The dis- 
persions in the inferred model parameters are a measure 
of the la random errors on the model parameters. In 
the Monte-Carlo simulations we also propagate the un- 
certainties on the model parameters that are kept fixed. 
We use A(to — M) = 0.10, Awgys = 0.6 km s~^ , and an 
uncertainty of 0.2° per coordinate in the plane of the sky 
on the SMC center position. 

LOS velocity studies provide some insight into the im- 
portance of rotation in the SMC. The old stellar pop- 
ulation of the SMC, traced by red giants, is consistent 
with a pressure-supported (V < a) spheroidal system 



with little evidence fo r rotation {Vrot < 17 km s~^ ; 
[Harris fc Zaritsk-vll2006[ ). This is consistent with the fact 
that many studies have f ound a large LOS depth in the 
SMC (e.g., ICrowl et al.l (f200Th 'l. By contrast the H I 
gas in the SMC shows rotation with a n amplitude of 
^ 40 km s~^ (Stanimirovic et al.ll2004[ ). This presum- 
ably indicates that the H I gas resides in a more disk- 
like distribution than the old stars. Young stars (spec- 
tral types O, B, A) show evidence for a rotation veloc- 
ity gradient of similar magnitude as the H I gas, but 
surprisingly, wi th a different major axis position angle 
(jEvans fc Howa rth 2008.^^ 

The results of .Harris fc Zaritskvl (|2Q06[ ) suggest that 
it is reasonable to assume that the SMC has limited ro- 
tation in the plane of the sky. If we analyze the PM 
data with fixed Kot = 0, the inferred SMC COM PM is 
Mw.LA/c = -0.772 ± 0.033 mas yr"^ and iin,lmc = 
— 1.117 ± 0.043 mas yr~^ . This can be compared to 
the result of just taking the weighted average of the 5 
SMC PM data-points, which yields fiw.LMC = -0.754 ± 
0.013 mas yr~^ and ^n,lmc = — 1.133±0.016 mas yr~-^ . 
These results are statistically consistent, despite the fact 
that the weighted average doesn't account for the in- 
fluence of viewing perspective. In Table [5] we list the 
estimate for the SMC COM that is obtained after sub- 
traction of the contribution from viewing perspective for 
each field. Taking the weighted average of results for in- 
dividual fields yields artificially low error bars, because 
it doesn't account for the actual scatter between mea- 
surements for different fields, or for uncertainties in the 
geometry and rotation of the SMC. 

To assess the importance of rotation, we have also per- 
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formed model fits in which the rotation is treated as a free 
parameter to be determined from the PM measurements. 
This yields K-ot = 29.5 ± 24.3 km s^^ . The PM data are 
therefore reasonably consistent with the absence of ro- 
tation, but they don't provide a very useful constraint. 
The corresponding best-fit SMC COM PM estimate is 
A*w,LMC = -0.694 ± 0.074 mas yr"^ and hn^lmc = 
— 1.055 ± 0.078 mas yr~^ . This is consistent with the 
estimate obtained assuming V^ot = km s~^ . However, 
the uncertainties are now significantly increased, due to 
the allowed variations in the rotation model correspond- 
ing to an amplitude variation AVJ-ot = 24.3 km . 

As our final SMC COM PM estimate we adopt the re- 
sult for Vrot ~ km s~^ , but allowing for a Icr Gaussian 
uncertainty AV^-ot- Based on the preceeding discussion, 
the average star used in our PM analysis probably rotates 
significantly less than the H I. We therefore adopt some- 
what arbitrarily AVrot = 15 km s~^ . This yields for the 
SMC COM PM that fiw,LMC = -0.772±0.063 mas yr"! 
and UN, LMC — — 1.117±0.061 mas yr~^ . The uncertain- 
ties are almost twice those quoted above, obtained with 
AK-ot = km s~^ . So as for the LMC, the quality of the 
data is now such that the uncertainties in the structure 
of the Cloud play an important role in how well we can 
estabhsh the COM PM. 

The value of (Xmin/^DOi=^)^/^ for the best-fit SMC 
model is ~ 2.8, independent of whether Vrot is treated 
as a free parameter or not. By contrast, for the best-fit 
LMC model derived in Paper II, {xl.iJNDOF)^^'^ = 1-8. 



Figure 7. (Top) Measured QSO motions (colored symbols-) with 
respect to stellar motions (open circles), similar to Figure [G] but 
now for the SMC. Green filled circles show results from all three 
epochs of data, and the pink squares for el and e2 only. The 
bottom panel assumes Vrot = km s~^ . 

Therefore, for both galaxies the measurements show 
somewhat more scatter than what can be accounted for 
by the models. This is not surprising, given the com- 
plexity of the data analysis, and the relatively ideal- 
ized nature of the models. Additional scatter might be 
due, e.g., to subtle detector geometric distortion varia- 
tions over time, residual astrometric effects due to imper- 
fect CTE, or more complicated galaxy structures than 
are accounted for in the models. It is intriguing that 
(Xmin/^ooF)^/^ is larger for the SMC than for the LMC. 
This might indicate that the structure and kinematics of 
the SMC are more complex (or less in dynamical equi- 
librium) than that of the LMC. 

4. COMPARISON TO PREVIOUS PM RESULTS 

4.1. Previous HST Results 

PM measurements from various sources, including 
the present paper, are summarized in Table |4l For 
the LMC, Kl obtained from two-epochs of data that 
Mm/.lmc = -2.03 ± 0.08 mas yr"^ and hn,lmc = 
0.44 ± 0.05 mas yr~^ . P08 obtained from the same data 
that ^J.w,LMC — — 1.956±0.036 mas yr^^ and ^n.lmc — 
0.435 ± 0.036 mas yr~^ . If we analyze the new PM data 
with the same fixed LMC center, inclination angle, and 
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line-of-nodes position angle from vdM02 used by Kl and 
P08, we obtain fiw.LMC — — 1.899 ± 0.017 mas yr^^ and 
fJ'N.LMC = 0.416 ± 0.017 mas yr~^ . This shows that the 
addition of our new third epoch of data, while keeping the 
LMC model the same, has the following effects: (a) the 
random errors in the final LMC COM PM components 
decrease, as expected based on the increased time base- 
line for about half the fields; (b) the value of hw,lmc 
decreases by about 1.6a compared to the earlier esti- 
mates; and (c) the value of fJ-N.LMC remains unchanged 
compared to the earlier estimates, within the uncertain- 
ties. 

Despite this good consistency, our final LMC result of 
MVK.LMC = —1.910 ± 0.020 mas yr^^ and fiN,LMC — 
0.229 ± 0.047 mas yr^^ deviates significantly from the 
final Kl and P08 estimates. This is because of the differ- 
ent LMC model used here, and in particular, the different 
position for the center. The difference is primarily evi- 
dent as a decrease by ~ 0.20 mas yr~^ in ijln,lmc- By 
contrast, our new value for iiw,LMC agrees with that of 
P08 at the l.lcr level. It is also worth noting that our fi- 
nal best-fit estimate does not have uncertainties that are 
much lower than those of P08, in particular in ^in,lmc- 
However, the reason for this is that P08 and Kl underes- 
timated the uncertainties, by not propagating all relevant 
uncertainties in the LMC disk model. Therefore, our new 
values are in fact significantly more accurate (both in a 
random and systematic sense) than the older ones. 

For the SMC, our results are more similar to those 
of P08 than K2. This is due primarily to the way in 
which we combined the results for different fields in K2. 
On a field-by-field basis, K2 and P08 actually obtained 
rather similar results. In K2 we didn't explicitly cor- 
rect for the astrometric effects of CTE. So we reported a 
value that was essentially an average of only 2 fields, be- 
cause 3 of the 4 used QSO fields had been taken with the 
same roll angle at both epochs. We also increased our 
errors to reflect a 'true' systematic plus random error, 
estimated from our analysis of the more numerous LMC 
fields. The K2 result therefore had large errorbars. While 
our new three-epoch result agrees with K2 for ij,n,smc, 
for fiw.SMC it differs by 2.4tT. Again, this is mainly be- 
cause of the way in which fields were combined in K2. 
The new three-epoch result, however, includes CTE cor- 
rections and a longer time baseline, and is therefore the 
more accurate one. 

Our three-epoch result for the SMC agrees reasonably 
well with that of P08. The value for hw,Smc agree within 
the uncertainties, while for fiN.SMC the results differ by 
1.7tT. This agreement is similar as for the LMC, if one 
compares analyses with the same assumed LMC center. 
So overall, there is acceptable internal consistency be- 
tween studies with different analysis methods, scientific 
instruments, and time baselines. 

4.2. Ground-based Results 

A number of older ground-based (and Hipparcos) anal- 
yses of the LMC and SMC PMs were discussed in Kl and 
K2. These studies had large uncertainties, but within 
these uncertainties the results were generally consistent 
with what we know now from HST data. 

There are now more recent results from two groups that 
are worth discussing (see Table |4]). iVieira et al.. (,201ft) 



measured the PMs of the Clouds based on CCD and plate 
material from the Yale/San Juan Southern Proper Mo- 
tion program, spanning a baseline of 40 years, and cov- 
ering a large area in the inter-Cloud region as well. They 
ultimately tie their reference frame to the ICRS. Within 
the uncertainties of ^ 0.3 mas yr~^ per coordinate, their 
LMC and SMC PM determinations agree with those pre- 
sented here. They obtain a stronger constraint on the 
relative motion between the Clouds, than on the mo- 
tion of either cloud individually: (/iiy , /XAr)sMC-LMC = 
(0.91,-1.49) ± (0.16,0.15) mas yr^^ . The correspond- 
ing difference between the Clouds' motions from our new 
HST results is: {i^w, I^n) sue -lmc = (1.188, -1.383) ± 
(0.039, 0.064) mas yr~^ . These results agree for ^jv, but 
differ by 1.7(t for jiw- Given the totally different method- 
ologies, this is a satisfactory result. However, the uncer- 
tainties on the ground-based data are large enough that 
they don't really help to validate the HST data at the 
level of its own uncertainties. 

Costa et al. performed two recent studies using the 
2.5m telescope at Las Campana s and centered on back- 
ground QSOs. In ICosta et al.l (2009) they measure 1 
fairly outlying QSO in the LMC over a 5 year baseline. 
The inferred LMC COM PM depends significantly on 
the assumed LMC rotation velocity. In Table 2] we quote 
the average of their results for Vrot = 50 km s"""^ and 
120 km s~^ (which corresponds to a rotation velocity 
consistent with that derived in Olseii ct al. (2 011 ) and 
Paper II). The uncertainties are 0.13-0.15 mas yr~^ per 
coordinate. Their PM value differs from our new HST 
result by 1.5a in fiw,LMC and 1.8cr in ijln,lmc- This is 
more than expected based on random errors alone. How- 
ever, we are unable to evaluate our HST results on the 
basis of this study, given that we have many more quasars 
a nd typically higher accuracy. 

ICosta et all ()201lD present an SMC COM PM estimate 
from 5 QSO fields over a 7 year baseline. This is similar 
to our study, but HST of course has better spatial reso- 
lution. Their PM value differs from our new HST result 
by l.ba in ^w,lmc, but agrees in ^n,lmc to within the 
uncertainties. Given the totally different methodologies 
compared to our work, this too is a satisfactory result. 
However, again the uncertainties on the ground-based 
data are large enough that they don't really help to val- 
idate the HST data at the level of its own uncertainties. 

5. GALACTOCENTRIC VELOCITIES 

The methodology for transforming an observed PM 
to a space motion in the Galactocentric rest frame 
is described in vdM02. The correction for the so- 
lar reflex motion requires knowledge of the Solar mo- 
tion in the Milky Way. In Kl and K2 we adop ted 
the standard lAU value (|Kerr fc Lvnden-Bel]iri986 i for 
the circular velocity of the Local Standard of Rest 
(LSR), Vq — 220 km s ~^ . However, models based 
on the PM of Sgr A* (jReid fc BrunthaleH 120041) and 
mase rs in high-mass star- formation regions ( Reid et al.l 
'2009^) have suggested that the circular velocity may be 
higher. iShatto w fc Loeb. (2009.) argued that this may 
signifi cant l y affe ct the orbit of the Magellanic Clouds. 
iMcMillanI (j2011f ) has presented a Milky Way analysis 
that includes all relevant observational constraints, from 
which he derived Vq = 239 ± 5 km s~^ . This is the 
value that we will adopt for the present study. We 
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also use the iMcMillaiil (|2011[ ) value for the distance of 
the Sun from the Milky Way, Rq = 8.29 ± 0.16 kpc. 
For the peculiar velocity of the Sun wit h respect to the 
LSR w e adopt the recent estimate from Sch onrich et al.l 
(M): (f/pec, VcW^pec) = (11.1,12.24, 7.25), with un- 
certainties of (1.23, 2.05, 0.62) km s~^ (being the quadra- 
ture sum of the random and systematic errors). Pre- 
vious work adopted th e solar peculiar velocity from 
iDehnen fc BinnevI (|1998[ ). However, there is now increas- 
ing evidence that Vpoc from that study is too small by 
^ 7 km s~^ . As a result, the solar velocity in the Galac- 
tocentric F-direction, vy = Vq + Vpcc, is 251.2 km s~-^ 
in our calculations here (consistent with the recent 
iBovv et al.l ()2012[ ) study who find vy 2421 J" km s'^ ). 
By contrast, it was 26 km s~^ lower in the calculations 
of Kl and K2. This directly impacts the LMC and SMC 
space velocity. 

In Table [S] we list the space velocity in the Galacto- 
centric rest frame of both the LMC (lines 1-6) and SMC 
(lines 7-10), implied by the various PM measurements 
available from HST. We also list the corresponding rela- 
tive velocity of the SMC with respect to the LMC (lines 
11-14). The relative velocity does not depend on the 
assumed solar motion. However, the individual LMC 
and SMC velocities do depend on the solar motion. For 
all PM measurements we list (labeled "lAU") the re- 
sults obtained under the assu mption of the lAU value 
of Vo = 220 km a'^ and the IDehnen fc BinnevI ()1998l ) 
solar peculiar velocity. For the new three-epoch mea- 
surements from the present paper we list also the results 
(labeled "new") obtained under the assumption of the 
McMill an (20ll value of Vq = 239 ± 5 km s~i and the 
Schonrich et all (|2010D solar peculiar velocity. For all re- 
sults we assume a distance uncertainty Am — M = 0.1 
mag. For the lAU value of Vq we also (arbitrarily) set 
the uncertainty to be AVq — 5 km s^^ ; the actual un- 
certainty must be significantly larger, if the McMillan 
(2011) results are indeed correct. The space motions 
that we derive based on the Kl, K2 and P08 PMs, with 
the "lAU" solar motion, are consistent with the values 
quoted in those papers. However, our errorbars are gen- 
erally larger and more accurate. This is because these 
previous studies did not propagate the error in distance, 
which is in fact a dominant error term (since v cx Dfi). 

The final results from the present paper are presented 
in lines 1, 7, and 11 of Table [5l This uses the three-epoch 
PM dataset. For the LMC, it uses the LMC geometry 
and rotation determined from the PMs in Paper II. The 
results are characterized by the following properties. The 
LMC has a total velocity (i.e., length of the Galactocen- 
tric velocity vector) of ?^tot,LMC — 321 ± 24 km s^^ . The 
SMC has a total velocity wtot,SMC = 217±26 km s~^. For 
both galaxies, the radial velocity z;rad is much less than 
the tangential velocity wtan- The value of Wrad,LMC > 0, 
so that the Clouds are past pericenter and moving to- 
wards apocenter. The relative velocity of the SMC with 
respect to the LMC is Utot,ici = 128 ± 32 km s^^ . The 
radial component of the relative velocity is almost twice 
the tangential component of the relative velocity. This 
implies a rather elliptical orbit, but a purely radial orbit 
(I'tan.rci — 0) is rulcd out. The value of Urad,rci > 0, so 
that the SMC is past its pericenter with respect to the 



LMC, and moving towards apocenter. 

The motions derived by Kl, K2, and P08 have been 
used for various studies in recent years. It is therefore of 
interest to examine how the new results differ from the 
old results, and why. 

Influence of solar velocity: Comparison of lines 2 and 
1 in Table [5] for the LMC, and lines 8 and 7 for the SMC, 
shows that going from the old to the new values for the 
solar velocity decreases all of Utot, Wiad and wtan, by 15- 
22 km s~^ . This is con sistent with arguments made by 
iShattow fc Loebl (|2009( ) about the influence of the solar 
velocity on the computation of the Clouds' orbits. 

Influence of LMC center position: Comparison of 
lines 3 and 1 shows that going from the vdM02 center 
(consistent with the center of the LMC bar) to the new 
LMC center from Paper II (consistent with the H I dy- 
namical center), reduces Wtot, Viad and Vt^n only mini- 
mally (~ 7 km s~^ ). However, ^j.n,lmc decreases by 
~ 0.2 mas yr~^ (see Table|4]) which equals ~ 47 km s~^ . 
This therefore changes primarily the angle of the tan- 
gential velocity as seen in projection from the Galactic 
Center, as previously pointed out by B07. The change is 
such so as to better align the past orbit with the location 
of the Magellanic Stream, but not to fully align it as we 
will discuss in fj6l 

Influence of third-epoch data for the LMC: Com- 
parison of lines 6, 5 and 4 shows that (when using the 
same solar velocity and the same parameters for the LMC 
geometry from vdM02) the addition of the third epoch 
WFC3 data and the reanalysis of the two-epoch ACS 
data decreases utot and Utan for the LMC by ~ 30 km s~^, 
when going chronologically from Kl, to P08, to the 
present paper. These changes are due to a progressive 
decrease in \^w,lmc\- The value of u^ad stays more-or- 
less the same. 

Influence of third-epoch data for the SMC: Com- 
parison of lines 10, 9 and 8 shows that (when using 
the same solar velocity) the addition of the third epoch 
WFC3 data and the reanalysis of the two-epoch ACS 
data decreases Vtot and t^tan for the SMC by ~ 65 km 
(a ~ la change from the K2 velocities and errors), when 
going chronologically from K2, to P08, to the present 
paper. These changes are due to how the K2 data from 
different fields were combined, which leads to a large dif- 
ference in the inferred nw component compared to what 
is obtained from a simple average of all the fields. The K2 
combination strategy also increased the error bars, how- 
ever, which is why there is consistency at the Itr level. 
The new data also indicate a small decrease in \fJ-N,SMc\ 
compared to P08. The value of Uiad changes by less than 
~ 15 km s^^ . 

Influence of third-epoch data on velocity uncer- 
tainties: Comparison of lines 6, 5, and 4, or 10, 9, and 
8, respectively, does not show much of a decrease in ran- 
dom errors on the velocities, when going chronologically 
from K1/K2, to P08, to the present paper. This is be- 
cause the PMs are now accurate enough that much of 
the velocity uncertainty is driven by distance uncertain- 
ties. These distance uncertainties were ignored in these 
previous papers, producing overly optimistic error bars. 

Many of these individual infiuences change the Galac- 
tocentric velocities in a similar direction. As a result, our 
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new velocities are quite different from those found by Kl 
& K2. 

New LMC velocity compared to Kl: Our new ve- 
locity has lower vtot and vtun by ^ 57 km s^^ . This is due 
to a combination of two effects: a ~ 30 km decrease 
due to a lower value of l/iw.LMcl and a ^ 27 km s~^ 
decrease to a revised understanding of the solar velocity. 
The latter also decreases Urad by ~ 25 km s~^ . More- 
over, we also derive a lower value of \hn,lmc\ due to 
a change in the derived LMC center. The new utot is 
approximately halfway between the traditional orbit re- 
quired by tidal Magellanic Stream m odels in a logarith- 
mic halo ([Gardiner fc Noguchil Il996l hereafter GN96), 
and the Kl result that led B07 to suggest using ACDM- 
motivated halo models that the Clouds may be on their 
first passage about the MW. We therefore perform new 
orbit calculations in SJHI 

New SMC velocity compared to K2: Our new ve- 
locity has lower Wtot and Wtan by ^ 85 km s^^ . This is due 
to a combination of two effects: a ~ 65 km s""'^ decrease 
due to a lower value of |/iH',s'Mc| (mainly due to how the 
fields were combined in K2) and a ~ 20 km s^^ decrease 
to a revised understanding of the solar velocity. The lat- 
ter is also primarily responsible for decreasing Uj-ad by 
^ 34 km s~^ . There is also a smaller contribution from 
the different SMC center used here compared to K2. 
New relative SMC-LMC velocity compared to 
K2: The total relative velocity wtot has changed very 
little. However, Uj-ad is now larger than Utan, so the orbit 
is more elliptical than before. The value of utot is not 
quite as high as found by P08, but consistent within the 
errors. 

6. ORBIT IMPLICATIONS 

The PMs and Galactocentric velocities from Kl & K2 
have been used in many studies to look at the past orbit 
of the Magellanic Clouds with respect to the MW, and 
the formation of the Magellanic Stream. As discussed in 
the preceeding sections, the latest HST observations and 
data analysis, combined with the latest understanding of 
the LMC geometry and solar velocity, imply decreased 
Galactocentric velocities, and improved understanding of 
the observational errors. We therefore rederive here the 
past orbit of the Magellanic Clouds based on the new 
results. 

6.1. Methodology and Parameter Space 

We explore a total of 5 models for the MW: 3 static 
models with total virial masses of 10^^, 1.5 x 10^^, and 
2 X 10^'^ Mq; and two cases in which we allow for the mass 
evolution of the MW as is expected from the hierarchical 
build up of mass in ACDM. While it is beyond the scope 
of this work to provide an exhaustive compilation of the 
observational determinations of MW mass, we do want 
to motivate our adopted range by pointing to the work 
of iGnedin et al.l (|2010( ). who obtain a MW virial mass 
of 1.6 X 1O^^M0 if the middle value for their data-based 
circular velocity determination is used. The uncertainty 
around this value is 20%. 

For the static MW cases we follow B07 and model the 
MW as an axisymmetric three-comp onent potential with 
a NF W halo, Miyamoto-Nagai disk fMivamo to fc Nagail 
|1975() and a Hernquist bulge (,Hernauist. ,19901^ In all 



cases the bulge is modeled with a scale radius of 0.7 
kpc and total mass of IO^^Mq. The NFW halo is adi- 
abatically contracted to accou nt for the p resence of the 
disk using the CONTRA code (|Gnedin et al. 2004). The 
NFW density profile is also truncated at the virial ra- 
dius (unlike in B07). While it is widely agreed that 
halos should be truncated, exactly how and where the 
truncation is done is arbitrary. The exact nature of this 
truncation does not affect our conclusions, because we 
consider a one-parameter family of models in which MW 
halo mass is considered a free parameter. 

Halo and disk parameters are listed in Table [6l In all 
cases the disk scale radius is kept fixed at 3.5 kpc while 
the total disk mass is allowed to vary in order to repro- 
duce the correct circular velocity at the solar circle. The 
MW circular velocity is taken as the update d value of 
239 k m s"i at the solar radius of 8.29 kpc (jMcMillanI 
120 lit ), rather than the standard lAU value, in all mod- 
els, unlike B07. Rotation curves for each MW model are 
presented in Figure [H We follow the methodology out- 
lined in B07 to capture the effects of dynamical friction 
acting on the Clouds owing to their passage through the 
MWs dark matter halo. 

In the cases where the LMC is not on a first pas- 
sage but has had a past pericentric passage, the typ- 
ical periods are quite large (> 5 Gyr). As such, we 
also want to investigate a simple model that accounts 
for the MW's mass evolution over the past 10 Gyr. In 
this case we model the MW only as a NFW halo (ne- 
glecting the disk, etc.). We model the mass evolution 
following the mean ni ass growth rate with redshift from 
iFakhouri et al.l ()2010D . The expected la scatter in this 
ev olution at each redshift is ro ughly 20% (from Figure 1 
of iBovlan-Kolchin et al. l '20101. We model the concen- 
tration evolution following Klypin et al. (2011), and the 
scaling of the viria l radius with redshift an d mass fol- 
lowing relations in iMaller fc Bullockl (|2004[ n. We con- 
sider current-day MW virial masses of 1.5 x lO-'^^M© and 
2 X 1O^^M0 in this exercise. We only explore the higher 
MW mass range because we are interested in constrain- 
ing whether the LMC is on a first infall (and as we will 
show below, it is always on a first infall in the low MW 
mass case). 

Our expectation is that the LMC mass is the dominant 
uncertainty in its orbital history (over MW mass evo- 
lution), since dynamical friction will change the LMC's 
orbit on timescales shorter than the MWs mass evolu- 
tiorQ. Therefore, for each static MW mass model, a va- 
riety of LMC masses are explored, ranging from 3 x 10^° 
to 2.5 X 10^^ Mq (the motivation for which is discussed 
directly below). This approach is different from B07, 
where we only considered one LMC mass of 3 x 10^'-' Mq. 

Our knowledge of the mass of the LMC is limited by the 
fact that kinematic data are only available in the inner 

^ The Fakhouri et al.l 1)20101 ') work is based on the Millenium 
II simulations and therefore uses an old value for ag , while the 
IKlvpin et all dMnTl work uses the Bolshoi simulations which have 
an updated value for ag. However, we do not think that this in- 
consistency should affect MW mass objects very much, since the 
old value of erg primarily suppresses power on low-mass scales. 

* However, the MW's mass evolution can substantially mod- 
ify solutions with o rbital periods of order > 6 Gyr, as the 
ll'akhouri et al.l 120101 ) relations imply t hat the MW's halo should 
be ~65% as massive 6 Gyr ago (see also lMoster et al]^012, ). 
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Figure 8. Rotation curves for the tiiree static MW potentials that we explore, with MW virial mass = 1 X (left), 1.5 X IO^'^Mq 

(middle), and 2 X (right). Vq, the Sun's velocity at the solar circle, is 239 km s~^ in all three cases, as indicated indicated by the 

horizontal dashed line in each panel. Rq, the solar radius, is kept fixed at 8.29 kpc, as indicated by the vertical dashed lines. The panels 
also indicate the amount of mass in the MW disk, as well as its disk scale-length. The horizontal black line shows the mean third-epoch 
LMC velocity derived here. The horizontal blue line shows the value for the MW escape velocity at 50 kpc. In the 2 X 10^^ Mq MW case, 
the escape velocity = 435 km s"'^ (just outside plotting range). 

~ 9 kpc. The observational estimat e of LMC total mass 
withi n this distance is 1.3 x lO^^M© (|van der Marel et al.l 
120091) . For each mass case considered here, the LMC is 
modeled as a Plummer sphere, where the softening radius 
is chosen such that the total mass contained within 9 kpc 
is ^ 1.3 X 10 ^ ^ Mp, a s observed, and as listed in Table [Sj 

iSaha et all ()2010f ) have detected stars out to 10 disk 
scale-lengths in the LMC (at a scale- length of 1.4 kpc, 
this amounts to stars as far out as 15 kpc). For a bound 
stellar component to exist at a minimum distance of 15 
kpc from the LMC center-of-mass, the tidal radius must 
be at least this large. In a 2 x IO^^Mq MW, the LMC's 
Roche radius at a MW-LMC separation distance of 49.5 
kpc reaches 15.4 kpc if the LMC total mass is SxIO^^Mq. 
This is the reason for our choice of lower bound in LMC 
mass. The highest total LMC mass explored is chosen to 
match the observed present-day LMC stellar mass and 
the re lations for the expected infall mass from lGuo et alj 
(|2010[) . For the evolving M W mass case we only consider 
LMC mass = 5 x 10^*^ M©. For this LMC mass we expect 
to find solutions where the LMC has completed an orbit 
about the MW, and this thus serves as a fiducial case for 
testing the plausibility of a non-first infall scenario. 

Whether the SMC is a binary companion to the LMC 
is primarily determined by the LMC mass, which pro- 
vides further motivation for varying the LMC mass. If 
the SMC is a long-term companion of the LMC, it will 
have been tidally truncated fairly early on, maintain- 
ing a roughly constant mass within its tidal radius at 
late times. Therefore, we do not vary the SMC mass. 
The SMC is modeled as a Plummer sphere of total mass 
3 X IO^Mq, with a Plummer softening of 1 kpc so that the 
mass within 3.5 kpc is 2.4 x IO^Mq (jStanimirovic et aLl 
120041 ). We do not account for dynamical friction acting 
on the SMC as it travels through the dark matter halo 
of the LMC. This is not significantly different than the 
approach used in B07. 

In our study we use the new 3rd epoch values and 
updated solar motion for both the LMC and SMC, and 
their corresponding errors (see lines 1 and 7 of Table [5]) • 
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Figure 9. The average number of pericentric passages that the 
LMC has made about the MW, as a function of LMC and MW 
mass, is indicated by different colors. Note that in all cases the 
LMC has at least made one pericentric passage since it is just at 
pericenter now, according to our data (N always starts at 1). The 
N=1.5 case represents orbits where there is roughly equal proba- 
bility that the LMC is on first infall (N=l) or has completed one 
orbit (N=2). The point size indicates the FIMS errors on the aver- 
age value, as indicated in the legend by black symbols. 



6.2. Orbits about the Milky Way 

As in B07 we follow the orbit of the LMC backward 
in time, using the present day velocities and positions 
and integrating the corresponding equations of motior0. 
For each combination of MW and LMC mass, we ex- 
plore 10,000 Monte-Carlo drawings from the LMC's ve- 
locity error distribution. We calculate the evolution of 
the LMC's Galactocentric radius, and keep track of how 

^ We implicitly assume here that the LMC and SMC are a binary, 
and that therefore the more massive LMC sets the orbital pat h. We 
ignore the SMC in this part of the calculation (but see ^6.311 . 
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many pericentric passages it makes about the MW over 
a Hubble time, assuming the MW mass is static over this 
time period. In Figure |9] we show the mean number of 
pericentric passages that the LMC has made about the 
MW, as a function of LMC and MW mass. Note that 
in ah cases the LMC has made at least one pericentric 
passage (N=l): our data indicate that the LMC is cur- 
rently just past pericenter, in agreement with previous 
works. This does not, however, imply that the LMC has 
made a complete orbit about the MW. The combination 
of LMC and MW masses that yield, on average, solu- 
tions where the LMC has completed at least one orbit 
are indicated by N> 2; first infall scenarios are indicated 
by N=l. N=1.5 represent cases where there is roughly 
equal probability that the LMC is on a first infall or that 
it makes one orbit about the MW (50% of the solutions 
have N=l and 50% have N=2). The size of the square 
indicates the RMS error on the mean value, as listed in 
the legend. 

For the majority of the cases we consider here, the 
LMC makes no additional pericentric passages about the 
MW than the one that it is currently at (it is on a first 
infall: dark blue squares). As might be expected, the 
LMC makes more pericentric passages if the MW mass 
is higher (larger binding energy) and the LMC mass is 
lower (minimal dynamical friction). The dispersion in 
the mean values increases as the mass of the MW in- 
creases; there are more solutions with one or more com- 
plete orbits allowed within the error space within a Hub- 
ble time for higher mass MW models. The dispersion de- 
creases as the LMC's mass increases; for high mass LMC 
mode ls (consistent with ACDM expectations, iGuo et all 
|2010| ) dynamical friction limits the number of closed orbit 
solutions within the error space. This forces the LMC on 
increasingly eccentric orbits with orbital periods longer 
than a Hubble time. 

Note that we have ignored the mass evolution of the 
LMC owing to the MW's tidal field in this analysis. Be- 
cause we have constrained all LMC models to match the 
total observed mass of the LMC within 9 kpc, higher 
mass LMC models will have more material at large radii 
and should thus be more affected by the omission of this 
physical effect than the lower mass LMC models. How- 
ever, Figure ^ illustrates that if the LMC today is well- 
described by the high mass models, it could not have 
made an earlier pericentric approach about the MW. As 
such, our argument is self-consistent. 

Figure |9] does illustrate that solutions can be found 
where the LMC makes more than one complete orbit 
about the MW within a Hubble time. For these cases. 
Figure [TOl illustrates the mean orbital period for the most 
recent orbit. Note that if the LMC has completed more 
than one orbit, the orbital period of the previous pas- 
sage is expected to be larger than the values quoted 
for the most recent passage because of dynamical fric- 
tion. For the 1 x IO^^Mq MW, the periods are > 9 
Gyr. For the 1.5 x lO^^M© MW, the typical periods are 
> 5 Gyr. In the most tightly bound scenario (LMC mass 
= 3 X lOi°M0, static MW mass ~ 2 x lO^^M©), the LMC 
completes on average 2 orbits about the MW where the 
most recent orbital period is 3-4 Gyr. This conclusion is 
no different than in B07 where we showed that the high 
MW model recaptures the isothermal sphere orbit from 
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Figure 10. The mean period is computed for only ttiose orbits 
that are closed within a Hubble Time, and is indicated using colors 
that correspond to duration. For example, if a particular MW- 
LMC mass combination produces 80% of LMC orbits with N=l, 
and 20% with N=2, then we only plot the average of the 20%. 
The symbol size indicates the RMS in distribution of periods, as 
indicated in the legend with black symbols. 

previous studies (albeit with larger orbital periods). 

For the cases where the LMC makes at least one com- 
plete orbit about the MW, i.e. LMC mass < 5 x lO^^M© 
and MW mass > 1.5 x lO^^M©, we look at the effect 
of a cosmologically motivated mass growth history for 
the MW on the orbital history of the LMC. In Fig- 
ure[TT]we show the LMC Galactocentric radius as a func- 
tion of time in the past {t = corresponds to today) 
for a current day MW mass of 1.5 x 10^'^ Mq (left) and 
2 X W^MQiright), and for LMC mass = 5 x 10^°Mq in 
both cases. In the first case, the period of the orbit is > 8 
Gyr for an evolving MW, which is ~ 2 Gyr longer than 
for a static MW. The LMC goes well outside the virial 
radius of the MW, i.e., the orbit is dramatically more ec- 
centric than in the static MW case. In the second case, 
it is shown that a 3 Gyr orbit would still on average in- 
crease by 1 Gyr if the mass evolution of the MW is taken 
into account. The period of the last complete orbit of 
the LMC about the MW is critical in figuring out the 
relative importance of MW tides in forming the Stream 
(discussed further in the next section). 

6.3. SMC-LMC Orbits and Implications for the 
Magellanic Stream 

We are also interested in re-evaluating the orbit of the 
Clouds about each other. A chance three-way encounter 
between the MW, LMC, and SMC at z = has low prob- 
ability, and therefore we seek to enumerate under what 
conditions, given the new velocities, the Clouds consti- 
tute a binary for a significant fraction of a Hubble time. 
Also, we have recently put forth a new Magellanic Stream 
model that relies on the past interactions between the 
Clouds when they are far from the MW virial radius, 
rather t han on the influence of the MW on th e Clouds- 
system (i Besla et al.ll2010l: iDiaz fc Bekkill2011[ ). The or- 
bital history of the SMC about the LMC is critical to 
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Figure 11. Left: the LMC Galactocentric radius as a function of time in tile past (t = corresponds to today) for a current day MW 
mass of 1.5 X IQI^Mq and an LMC mass of 5 X 10^" Mq. The period of the orbit is 2 Gyr longer for an evolving MW (solid line), and 
the LMC goes well outside the virial radius of the MW, i.e. the orbit is dramatically more eccentric than in a static MW (dashed line). 
Right: the same for a current day MW mass of 2 X 10^^ Mq. The shortest orbital periods we can obtain are ~ 4 Gyr. 

assess the viability of such a model. 

We draw 10,000 cases from the LMC and SMC PM 
error distribibutions in Monte-Carlo fashion, and orbits 
are computed for each combination of LMC and MW 
masses. The SMC is assumed to be tidally truncated 
and hence its mass is kept fixed at 3 x 1O^M0. At each 
time-step, we compute the escape velocity of the SMC 
from the LMC, Vcsc, and compare this to the relative 
velocity between the Clouds, Vls- We do not require 
them to be bound today, but rather search for cases at 
some point in time where Vls < Vcsc, and keep track 
of the longest amount of time over which this condition 
remains true. 

This analysis, therefore, accounts for capture events 
and the fact that the MW's tidal field might be disrupt- 
ing the binary today. By also requiring that the escape 
velocity condition be satisfied, as opposed to searching 
for minima in the distance between the Clouds as many 
previous authors have done, we are choosing real bi- 
nary configurations, rather than chance superpositions 
caused by cases in which the SMC is on its own ellip- 
tical orbit about the MW, and therefore formally gets 
closer and farther away from the LMC even though the 
LMC may be on a very different orbit (as illustrated in 
iRilzicka et al|[2010t ). 

In all searched cases the SMC has made at least one 
close encounter with the LMC in the past. This is unsur- 
prising, as it is a direct result of the relative orientation 
of their three-dimensional velocity vectors. The SMC is 
currently ^ 20 kpc from the LMC and moving away from 
it. Integrating the SMC's orbit backwards in time will 
thus always yield a close encounter with the LMC within 
the past 500 Myr. This conclu sion has bee n drawn by 
many previous authors (Ruzicka et al.llTolOl GN96) and 
is the basis of the LMC-SMC collision theory put forth 
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Figure 12. The percentage of cases from our 10,000 Monte Carlo 
drawings in which Vj^s < ^csc for the SMC from the LMC for more 
than 2 Gyr in the past, shown as a function of LMC and MW 
mass. The legend shows the color-scheme indicating the range of 
percentages for each combination of LMC and MW mass. 



in iBesla et al.l (goH). 

In Figure [T^] we show the outcome of our search for 
binary orbits between the Clouds. The percentage of 
Monte Carlo drawings for which Vls < Vesc for more 
than 2 Gyr in the past is shown as a function of LMC 
and MW mass. We have chosen the 2 Gyr duration since 
this is the estimated age of the Magellanic Stream (as 
discussed below), and so the LMC and SMC should be 
a binary for at least this amount of time in each LMC 
and MW mass combination. Colored squares indicate the 
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corresponding percentages, as indicated in the legend. 

Unsurprisingly, the percentage of cases for which Vls < 
Vcsc increases as the mass of the LMC increases (increas- 
ing binding energy). However, we also note a dependence 
on MW mass. At an LMC mass of 3 x 10^°Mq, the per- 
centage of orbits that satisfies this criterion is 20% if the 
MW mass is 1 x lO^^j^^^^ xhis number drops to 4% if 
the MW mass is 2 x 10^"^ Mq, implying that bound con- 
figurations between the LMC and SMC are statistically 
unlikely for low LMC mass and high MW mass. As the 
LMC mass increases, these statistics improve, but for a 
MW mass of 2 X IO^^Mq, they only reach 20% once we 
get to a LMC mass of 8 x 10^°Mq. The corresponding 
percentage for a I x lO^^M© MW is 50%. 

The RMS of the distribution of cases matching this 
escape velocity criteria increases as the LMC mass in- 
creases, representing a larger spread in the possible du- 
rations. As discussed above, as the MW mass increases 
there are much fewer binary LMC-SMC orbits, even for 
high LMC mass. The PMs indicate that the current rel- 
ative velocity between the Clouds is high. This implies 
that the SMC must be on an eccentric orbit about the 
LMC. For larger MW mass models, the tidal field is more 
effective at disrupting such a wide binary pair. 

Therefore, from our statistical analysis we conclude 
that large LMC masses 1 x 10^ ^Mq) are favored if 
the Clouds are to have been in a binary, and that fur- 
ther, a MW mass < 1.5 x 10 ^^M0 is also required. From 
Figure |9l orbits in this mass range almost always cor- 
respond to first infall scenarios. The assumption that 
the Magellanic Clouds constitute a long-lived binary pair 
thus implies that the Clouds are likely on their first infall 
about the MW. 

Based on our searched parameter space, and the re- 
quirement that the LMC and SMC have been a long-lived 
binary, we adopt a canonical model with the highest MW 
mass and lowest LMC mass that will fulfill these crite- 
ria. This gives (LMC, MW) mass = 1.8 x IO^Mq, 1.5 x 
IQi^Mp. 

In Figure [T3] we show an example of one of these long- 
lived binary states (N=4-6) in the canonical model (note 
again that the SMC is assumed to be tidally truncated 
and therefore kept at fixed mass in our models) . The left- 
hand panel shows the Galactocentric radius of the LMC, 
the Galactocentric radius of the SMC and the relative 
distance between the two Clouds as a function of time. 
The middle panel shows the past orbits of the Clouds in 
the Galactocentric {Y — Z)-plane. The right panel shows 
the past orbit as a function of {l,b). Interestingly, the 
new velocities now cause the past SMC orbit to cross 
over in the direction of the Stream, something that we 
did not achieve with the K2 velocities, but was advocated 
for in earlier works like GN96. 

We note that with the vdM02 center, the relative LMC- 
SMC velocity is 143 ± 31 km s^^ (Table O. With the 
new center adopted in this work, it is 128 ± 32 km . 
Since the relative velocity is higher in the former case, the 
orbit of the SMC would be less bound at fixed mass, and 
therefore we expected fewer bound cases. Suprisingly, we 
in fact find no binary LMC-SMC orbits when using the 
vdM02 velocities. We attribute this not just to the larger 
absolute value of relative velocity, but to the fact that the 
angle between the velocity vectors is larger. This is an 



interesting finding that argues in favor of the new center 
used here and derived in Paper II. 

While dynamical friction owing to the passage of the 
Clouds through the MW's dark matter halo was ac- 
counted for, we did not explicitly account for dynamical 
friction acting on the SMC owing to its passage through 
the LMC's dark matter halo. The inclusion of this effect 
would place the SMC on an increasingly eccentric orbit 
in the past. We have recently advocated for a high eccen- 
tricity SMC orbit about the LMC in our simulations of 
the Magellanic Stream; such eccentric orbits are needed 
to explain the large angular extent of the Stream an d 
to prevent the Clouds from merging (|Besla et al.llMol) . 
However, such eccentric orbits are more easily disrupted 
by the MW's tidal field and thus work against main- 
taining long-lived binary LMC-SMC states. This further 
strengthens our argument that the preferred mass scale 
for the MW is less than 2 x IO^^Mq, since in reality there 
should be even fewer viable binary orbits in this model 
than we currently find. 

Another way to approach these questions is to look 
in more detail at the LMC-SMC dynamics in LMC-MW 
mass combinations that give rise to more traditional or- 
bital trajectories for the Clouds about the MW. Specifi- 
cally, consider the case of a low-mass (tidally truncated) 
3 X 10i°A/o LMC orbiting around a 2.0 x lO^^M© MW, 
with an orbital period of ~ 4 Gyr. We draw 10,000 
LMC and SMC PMs in Monte-Carlo fashion from their 
error distributions, and plot the relative proper motion 
between the two Clouds in Figure [141 Each value is 
color-coded by the length of time for which the crite- 
rion Vls < Vesc held true in the past. While it is clear 
that the large majority of pairs within the la error el- 
lipse do not produce a bound system for > 2 Gyr, it is 
interesting that some orbits with a bound pair in excess 
of several Gyr do exist within this error-space. A larger 
abundance of such orbits can be found at relative PMs 
that are more in the West direction than our measure- 
ments imply. So the combination of galaxy masses shown 
in this plot, while unlikely, cannot be ruled out strictly 
based on the requirement that the LMC and SMC must 
have been bound for at least 2 Gyr. However, it should 
be noted that this does not imply that traditional models 
for the Magellanic Stream (which have often used masses 
such as shown in this figure) are tenable, because such 
models produce orbital periods that are inconsistent with 
the age of the Stream (discussed more below). 

In summary, perhaps unsurprisingly, we find that to 
obtain long-lived binary LMC-SMC orbits the LMC mass 
must be relatively high, and the MW mass must be rel- 
atively low. The lower mass of the MW implies that the 
Clouds are on a first infall and, as such, the MW's tidal 
field would be insufficient to disrupt the binary pair. We 
are more readily able to find such long-lived orbits with 
the new LMC center derived here rather than that pre- 
viously derived by vdM02. We also find that a close en- 
counter between the Clouds in the recent past is a generic 
result. 

Given that the Clouds are a binary pair, the more 
massive LMC sets the orbital period for the Clouds 
about the MW. As shown above, the LMC mass must be 
> 1 X IO^^Mq for the Clouds to have been a long-term 
binary system. The smallest orbital periods about the 
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Figure 13. Le/i; Evolution of the LMC (red) and SMC (green) Galactocentric radius as a function of time in the past, for our canonical 
orbit, i.e., Mvir = 1.5 X lO^^M©, Mlmc = 1-8 x lO^^M©, Msmc = 3 x 10^ Mq, and the mean third epoch LMC velocity. The blue line 
shows the relative distance between the two Clouds. The dot-dash horizontal line shows the virial radius, iJvin for this MW model. Middle: 
the same shown in the Galactocentric {Y, Z)-plane. The orbits of the LMC (solid red line) and SMC (dashed red line) are followed for the 
past 4 Gyr, and their present locations with respect to the MW are marked with black squares. The MW center is marked and its disk 
orientation is indicated by the short dash. iJvir is also shown. Right: the same orbits for the LMC (solid red line) and SMC (dashed red 
line) are shown in the (I, f))-plane. The blue lines show the corresponding GN96 orbits, which were chosen by those authors to match the 
location of the Hi in the Magellanic Stream (see Figure 9 of B07), and stands in for the Stream location here. 
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Figure 14. A further exploration of LMC-SMC orbits in a low- 
mass LMC and high-mass MW model (masses indicated in the 
figure). This combination of masses gives rise to more traditional 
LMC orbits (albeit with longer orbital periods as can be seen in 
Figure [TTJ. We did 10,000 Monte Carlo drawings from the PM 
error distributions. We show the relative PM between the Clouds 
with each point color-coded by the amount of time for which Vls < 
VcBc in the past (colored squares). The black triangle is the mean 
relative PM value from this study, with la, 2a, 3a, and 4a error 
ellipses indicated as well (dashed). While the vast majority of 
drawings do not give rise to orbits that were bound for more than 
2 Gyr (cyan squares), there are some orbits that are bound for 2-4 
Gyr within the la error ellipse (dark blue squares). 



MW that we can obtain for such large LMC masses are 
in excess of 4 Gyr. This is problematic for models that 
rely on Milky Way tides to form the Stream, because age 
estimates of the Stream imply that the Stream is a young 
feature, much younger than 4 Gyr old. In fact, given the 
rate at which the Stream is being ablated away, as in- 
ferred from simulations of the observed anomalously high 



Ha emission in the S tream f Weiner fc Williamsl[T996l ) by 
[Bland-Ha wthorn et a l. (2003), and based on estimates of 
the survivabilit y of high velocity clouds in MW-type envi- 
ronments ((Heitsch fc PutmanI 120091 : iKeres fc Hernguis^ 
I2009f ). it is unlikely that the Stream could have survived 
about the MW for more than ~l-2 Gyr. 

The past orbit of the LMC and SMC, implied by the 
new velocities, is now better aligned with the location 
of the Stream HI than it was with the K1/K2 veloci- 
ties. As discussed in B07, this is directly due to the 
lower value for hn, which in turn is obtained here di- 
rectly from the dynamical center of the LMC having 
changed with respect to what was used in Kl. How- 
ever, we note that the LMC orbit still does not trace 
the Stream closely: there remains a sizeable offset. The 
new SMC PMs, however, do allow the past SMC orbit to 
cross over the Stream location, which is now easier to rec- 
oncile with models where the Stream is largely stripped 
from the SMC (GN 96; Rtizicka ct al. 2010; Diaz & BekO 
120111: IBesla et al.|[2012i) than it was with the K2 SMC 
PMs. 



6.4. First Infall? 

As we learned from B07, the question of first infall 
is a model-dependent rather than simply a velocity- 
dependent one. We showed in that work that no mat- 
ter whose velocities were used (Kl/GN96/prior ground- 
based determinations), the Clouds were always on a first 
infall in a 1.2 x lO^^^fg NFW MW halo, and that in or- 
der to recapture something resembling the 'traditional' 
isothermal sphere orbit, we had to go to a 2 X lO^^Mg 
NFW halo. 

Given that the new Galactocentric velocities derived 
here are lower than those in Kl, we want to critically 
re-examine the arguments for whether the LMC may be 
on a first passage, which we define as an orbital solu- 
tion wherein the LMC has first entered the virial radius 
of the MW within the past 1-4 Gyr and has not com- 
pleted an orbit in that time. We have already shown 
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above that with our new velocities, a low LMC mass 
can make past pericentric passages about a high mass 
MW. Therefore, arguments for a first passage, being stiU 
model-dependent in nature, can be cast as arguments 
against a low-mass LMC and a high-mass MW. Our ar- 
guments fall into three main catagories: (1) orbital ec- 
centricity and cosmological expectations, (2) LMC tidal 
radius, and (3) LMC-SMC binarity. We also discuss a 
fourth argument that attempts to reconcile the high gas 
content and active star formation of the Clouds with their 
current close proximity to the MW. 

Orbital Eccentricity and Cosmological Expecta- 
tions: As shown in Figure[TU]in all cases where the LMC 
does complete an orbit about the MW, the corresponding 
orbital periods are typically large (> 4 Gyr). Such orbits 
take the LMC to apocentric distances of order the virial 
radius of the MW and thus imply large orbital eccentric- 
ities. The computed eccentricities for the 3 x IO^^Mq 
LMC and the 2 x 10^°Mq evolving MW model are 0.6- 
0.7. These eccentricities get hi gher for larger LMC 
masse s. Comparing to Figure 5 of iBovlan-Kolchin et al.l 
(|2011h . only 20% of LMC analogs accreted at early times 
(> 8 Gyr) have such orbits. Indeed, in the prevailing 
ACDM model of hierarchical structure formation, mas- 
sive satellites that are accreted at early times are very 
rarely found o n highly eccentric orbits about MW-type 
hosts at z=0 (IBovlan-Kolchin et al.|[20Tlt iWetzell [2010l: 
iStewart et al.l l2008f). 

Generally, subhalos that are accreted on ^radial or- 
bits at early times are either preferentially d estroyed 
or ex ist on more circularized orbits today (e.g., lBenson| 
120051 ). While 50% of LMC analogs accreted more than 
4 Gyr ago and 60% of LMC analogs accreted more 
than 8 Gyr ago have eccentricites < 0.5 in the Mil- 
lenium 11 sample, only 20% of LMC analogs accreted 
within the last 2-4 Gyr hav e such modest eccentricites 
(Boylan-Kolchin et all 1201 ID . A high MW mass, low 
LMC mass today isn't ruled out from cosmological expec- 
tations, but the high eccentricities (0.6-0.7) calculated 
from our orbits are more typical of LMC analogs accreted 
within the past 4 Gyr. Therefore, from the eccentricities 
of our orbits, a first infall scenario for the Clouds is the 
favored orbital solution of ACDM theory. 

We can also make a more general timing argument 
based on the fact that the Clouds and the MW must 
have been in close proximity at the time of the Big Bang 
(i n analogy with similar arguments for e.g., M31 and Leo 
I (|Li fc W hite 2008)). This implies that the Clouds had a 
pericentric approach with the MW at the time of the Big 
Bang. So if there has been more than one complete orbit 
since, then the period must be < To/2, where Tq = 13.73 
Gyr. Hence, orbits with a pericenter > 6.9 Gyr ago are 
not physical, because the Clouds and MW were not to- 
gether at the Big Bang. Of course, this is oversimplified 
for many reasons, e.g., the mass of the MW increases 
with time. But as we show in Figure 1111 this typically 
increases the period by ~ 2 Gyr. Therefore, finding an 
orbit that is consistent with the proper motion data with 
a pericenter at e.g., 10 Gyr ago, does not rule out a first 
infall scenario, since such orbits are not physical based 
on such a timing argument. 

LMC Tidal Radius: As discussed in ^Othc L MC to- 
tal mass must be at least 3 x lO^^M© from the .Saha et al.l 



(|2010f ) observations. iMunoz et al.l ()2006( ) have claimed 
a detection of LMC stars out to 20 kpc. If stars do 
indeed exist out to 20 kpc, the required LMC mass is 
^ 8 X 1O^°M0. It should be pointed out th at more tradi- 
tional Magellanic Stream models (GN96; iDiaz fc Bekkil 
,2012 ) assume a total LMC mass of 1 x 10^°Mq which 
is already ruled out by the observations. Their choice 
of LMC mass is important to point out since this is the 
reason their orbital periods are so different from ours. 

LMC-SMC Binarity: Studies of the star-formation 
histories of the Clouds all agree on the fact that the 
star-formation rate started to increase ~ 4-6 Gyr ago, 
and was q uiescent in bot h galaxies before this time 
ijHarris fc Zaritsky 2009 : Cignoni et al.ll2012D . The co- 
incident increase in star-formation in both galaxies may 
imply that they were in a common envelope at the time, 
and correspond s to the epoch th at the LMC captured 
the SMC in the lBesla etaH (j2012h model. 

If we therefore take the stance that the two Clouds 
must have fallen into the MW as a bound system (al- 
though they need not necessarily be bound at present) 
we find that we are able to place limits on both the MW 
and LMC galaxy masses. We find that in general an 
LMC mass > 1 x lO^M© is needed for the SMC to have 
been bound to the LMC for > 2 Gyr in the past, and 
that further a MW mass < 1.5 x lO^^M© is required. 
Note that LMC orbits at this particular mass are almost 
always on first infall regardless of MW mass. The LMC 
mass required to keep the SMC bound to it goes up as the 
mass of the MW increases, in order to compensate for the 
increasing tidal field of the MW. But as the LMC mass 
increases the orbital eccentricity does as well, making a 
first infall scenario more likely based on the plausibility 
argument above that high eccentricities and early infall 
times are relatively rare in cosmological simulations. 

It is harder to maintain a LMC-SMC binary for long 
periods of time in high mass MW and low mass LMC 
models: less than 5% of the 10,000 searched Monte Carlo 
cases result in binary configurations, meaning such mod- 
els generally require that the LMC capture the SMC 
while in orbit about the MW - a statistically improb- 
able event. This argument implies that the MW's mass 
is likely less than 2 x IO^^Mq. 

The Clouds and the M orphology-Distance rela- 
tion: Ivan den Berghl ()2006r ) argued for a recent accretion 
of the Clouds, based on a comparison of the morphologies 
of MW and M31 satellites that showed that the Clouds 
are the only gas-rich dirr galaxies at small Galactocentric 
distances, all other dirrs being at large distances from 
their respective hosts. Some recent studies have looked 
at this issue from a larger statistical/cosmological point- 
o f-view. 

iToUerud et al.l (|2011h study a volume-limited spectro- 
scopic sample of isolated galaxies in SDSS and find that 
bright satellite galaxies around MW-type hosts are sig- 
nificantly redder than typical galaxies in a similar lumi- 
nosity range, and argue that this is indicative of environ- 
mental quenching. This is found to be in stark contrast 
to the LMC, which is anomalously blue in comparison to 
other LMC-MW analogs in SDSS. The authors attribute 
this to the fact that the LMC may be undergoing a trig- 
g ered star f ormation event upon first infall. 

iGeha et al.. (2012,) have used the NASA-Sloan Atlas 
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to demonstrate that dwarf galaxies in the field (with 
masses in the LMC range) all have active star-formation 
(< 0.06% of field-dwarfs in their sample have no star- 
formation). By contrast, the majority of quenched galax- 
ies are all within 2-4 virial radii of a massive host. There- 
fore, ending star-formation in such dwarf galaxies ap- 
pears to require the presence of a more massive neighbor. 
Had the Clouds been accreted early in the Universe then, 
it is unlikely that these galaxies would currently have as 
much gas or sta r form ation as they do. 

Wetzel ct al. (|2012[ ) used SDSS galaxy group/cluster 
catalogs together with N-body simulations to look at the 
quenching timescales of satellites at z ~ 0. They find a 
quenching scenario in which satellite star-formation his- 
tories are unaffected for 2-4 Gyr after infall but then 
quench rapidly (with an e- folding time of < 0.8 Gyr). 
Interestingly, because of the time delay before quenching 
starts, satellites are found to experience significant stel- 
lar mass growth after infall, which the authors point to 
as a key reason for the success of the subhalo abundance 
matching technique: supporting the larger LMC masses 
that we favor in this study. 

In conclusion, despite the revised lower estimates for 
the 3D velocities of the Magellanic Clouds presented in 
this work, taken together. Figures El [101 HU and [T2] make 
a strong case for a first infall scenario. This conclusion 
draws largely from a recognition that the fundamental 
change to our understanding of the orbital history of the 
Clouds comes not only from the PMs, but also from a 
cosmologically motivated understanding of the mass evo- 
lution and dark matter halo profile of our MW. This pic- 
ture is thus consistent with the theory that the origin of 
the large scale gaseous structures of the Magellanic Sys- 
tem (Stream, Bridge and Leading Arm) and the internal 
structure and kinematics of the Clouds are a result of in- 
teractions between the LMC and SMC in a fir st infall sce- 
nario, rather than interactions with the MW ()Besla et al.l 
[2bl0 . 2012 ). 

7. CONCLUSIONS AND DISCUSSION 

We have analyzed a new third epoch of data of 10 
QSOs behind the LMC and 3 behind the SMC with 
WFC3/UVIS (one LMC field, L22, has only first epoch 
ACS and third epoch UVIS data). We combine these 
data with previously obtained ACS/HRC data giving a 
~ 7 yr baseline, in order to measure PMs for both Clouds, 
and to refine and validate measurements that used only 
two epochs over a 2 yr baseline (Kl, K2, P08). We have 
also reanalyzed the first two epochs of data for 21 LMC 
fields and 5 SMC fields, using a method to account for 
the degrading CTE of the HRC as well as an improved 
approach to linear transformations. Here we summarize 
our findings, the implications for the Clouds' orbits, and 
provide a look ahead to the future. 

7.1. Conclusions for COM PMs and Calactocentric 
Velocities 

The ACS reanalysis and the third epoch analysis give 
very consistent results at the level of the measured per- 
field PMs. In order to obtain COM PMs from these 
field PMs, the geometry and internal motions of the 
L/SMC need to be taken into account. The addition 
of the WFC3 epoch provides very small per-field PM er- 
rors of ~ 0.03 mas yr~^ (7 km s~^ , or 1.6% of the total 



PM). We are therefore able to independently constrain 
all parameters of the LMC PM rotation field, including 
the inclination, i, the position angle of the line-of-nodes, 
9, the dynamical center, rotation velocity amplitude, and 
even the distance. This is a big improvement over Kl. 
This procedure and the results are discussed in detail in 
Paper II. The features of import for the present study are 
that we obtain a rotation velocity t hat agrees ve r y weh 
with that obtained independently bv lOlsen et al.l (201 If ) 
using a very large number of LOS velocities of stellar 
tracers. Moreover, we find strong evidence that the cen- 
ter of the LMC (as obtained from the PM fit to the ro- 
tation field) is consistent with the HI dynamical center, 
and not the center derived by vdM02, which agrees with 
the brightest part of the LMC bar. 

For the SMC, the sparse coverage of QSO fields ren- 
ders it difficult to gain much insight into its geometry. 
We therefore fit a relatively simple model to the SMC 
PM data, which allows for viewing perspective and a sin- 
gle overall rotation of the SMC in the plane of the sky. 
The old stellar pop ulation in the SMC show s little evi- 
dence for rotation (Har ris &: Zaritskvl |2006( ) so we base 
our analysis of the SMC COM PM on the assumption 
that Vrot = ± 15 km s~^ . Analyses in which Vrot is 
instead fit to the PM data do not yield a significantly 
different answer. 

If we analyze the new LMC PM data with the same 
fixed geometric parameters as used in vdM02 we obtain a 
COM PM whose (i) value for fiw is lower, but consistent 
within 1.6a with the Kl result, (ii) value for fi^ remains 
unchanged with respect to Kl, (iii) random errors are 
smaller by a factor of almost 5 than Kl. If we analyze 
the PM data with the LMC model derived from the data 
itself, which is the preferred approach, then the value 
of fiw and its random error remain largely unchanged. 
However, the value for /i^r changes significantly (by 4ct), 
and its random error increases by a factor of almost 3. 
For the SMC, the three-epoch analysis agrees with K2 
for fiN, while for fiw it differs by 2Aa. This is largely 
due to how fields were combined in K2 and not due to 
intra-field differences. Our final SMC COM PM value is 
in rough agreement with the results of P08. 

We attempt here to give as accurate an estimate as 
possible of the observational random errors in the COM 
PM, by propagating all unknowns in the geometry as well 
as the PM determination. This obviously increases errors 
with respect to previous studies which only propagated 
PM errors. Remarkably, the PM data are now no longer 
the dominant source of uncertainty, but rather uncertain- 
ties in the structure of the Clouds, dominate how well the 
COM PM can be established. Our final COM measure- 
ment errors, listed in line 1 of Tabled do not represent a 
huge improvement over what was listed in previous HST 
works, despite the significantly improved accuracy of the 
present work. However, this is because those previous 
works underestimated the true random errors. Regard- 
less, the uncertainty in the inferred transverse velocities 
of the Clouds are now dominated by their distance un- 
certainties, and not their COM PM uncertainties. 

In order to turn the measured COM PMs into Calac- 
tocentric velocities, we need to know the solar veloc- 
ity. Our understanding of the solar velocity has re- 
cently been revised upwards, and this works to directly 
reduce the Calactocentric velocities of the Clouds by 
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a comparable amount. Therefore, the Galactocentric 
vtot = 321 ± 24 km s^^ for the LMC presented here is 
57 km s""'^ lower than that in Kl, owing roughly equally 
to the decrease in fiw described above and to the increase 
in solar velocity. There is a much smaller dependence 
on the new dynamical center derived here (see Table [5]). 
The new center affects mainly the /^at value, which de- 
termines the location of the orbit in projection on the 
sky rather than the tangential velocity (as described in 
B07). By contrast, if we use the same geometric model 
and solar velocity (lAU value) used in Kl, wtot decreases 
by 31 km compared to the Kl value, which is com- 
parable to the size of the la error bar. 

The SMC's inferred Vtot = 217 ± 26 km s^^ is 
85 km s~^ lower than in K2. This is due mainly to the 
new proper motion derived here (which itself is due to 
how the fields were combined in K2, rather than differ- 
ences in the derived per-field PMs). Roughly 20 km s~^ 
of the decrease in vtot is due to the revised solar velocity, 
and there is also a small contribution from the different 
SMC center used here, compared to the K2 analysis. 

7.2. Conclusions for the Clouds' Implied Orbits 

Given the new velocities obtained here, we have re- 
evaluated the past orbital histories of the Clouds about 
the MW. We find that the dominant unknowns are the 
MW and LMC masses. However, some reasonable ar- 
guments allow us to narrow down the allowed ranges of 
MW and LMC masses. The results continue to make a 
strong case, as first argued in B07, that the Clouds are 
likely on their first infall into the MW. 

From the search for bound orbits between the two 
Clouds (see Figure [T^ , we are able to identify a com- 
bination of masses that satisfy the criterion of long- 
lived binarity, and also produces orbits around the MW 
that are plausible from a cosmological point of view 
(jBovlan-Kolchin et al.l[201ll : iBusha et aLlboilal) . This 
yields a preferred MW mass, < 1.5 x 10^^ M0. The 
probability of a stable LMC-SMC binary configuration 
decreases as the MW mass increases. An LMC mass 
> 1 X 1O^^M0 is needed in order to keep the SMC 
bound to the LMC for a reasonable fraction of a Hubble 
time. Taken together, this combination of MW and LMC 
masses imply that the LMC / the Clouds are on their first 
infall (see Figure [S]). 

From a cosmological point of view, these orbits do not 
represent a major alteration to the interpretation given 
in B07. Here we have explored a more full set of mass 
models for the LMC. The inferred LMC mass agrees with 
the expected total infall mass of the LMC, as calculated 
from the observed present-day stellar mass (v dM02) and 
the r elations fro m ha lo occupation models (|Guo et al.l 
l2010t [Sales et al.ll2011h . A generic outcome of our study 
is that regardless of what MW mass is used, large LMC 
masses are always on a first infall. 

We have also implemented a simple model for the ex- 
pected mass evolution of the MW over the past ~ 10 Gyr. 
This, as might be expected, significantly increases the pe- 
riods of LMC orbits that do make a previous pericentric 
passage (higher MW mass models), i.e., the orbits are 
highly eccentric. In a recent study, motivated by the dis - 
cussion over the higher Solar velocity. iZhang et al.l (|2012[ ) 
have also revisited the Clouds' past orbital history, uti- 



lizing numerical simulations to explore the evolution of 
the MW. Like us here, they find that it is possible for the 
LMC to make past pericentric passages given the obser- 
vational and theoretical error-space. They do not vary 
LMC mass (it is kept fixed at 2 x 10^° Mq), or investi- 
gate whether binarity imposes constraints, so we cannot 
directly compare our results with theirs, but were we to 
use such a low LMC mass, our studies would likely be 
consistent. 

We put forth four arguments why we think orbits in 
which the LMC/SMC make a previous pericentric pas- 
sage are implausible. One is the expected eccentrici- 
ties of orbits from cosmological simulations along with 
a simple version of the timing argument (see § 16. 4[) , the 
second is the tidal radius of the LMC, and the third is 
the fact that long-lived SMC-LMC binary configurations 
have lower probability for such orbits, given the observed 
COM PMs of the Clouds. The fourth a rgument derives 
from the work of Ivan den Bergbl ()20Q6f ) who conducted 
a morphological comparison of the satellites of the MW 
and M31, finding that the L/SMC are the only two gas- 
rich dirrs at close Galactocentric distance to a host. He 
argues that they may be interlopers from a remote part 
of the Local Group rather than true satellites of the MW. 

This line of argument is also consistent with the fact 
that LMC-SMC-M W analogs are r elatively rare today in 
our local volume ([Liu et al. '20 1 11) and in cosm ological 
simulations at z=0 (Bovlan-Ko lchin et al.l 1201 It ). MW- 
type hosts are efficient at tid ally disrupting s uch b ound 
configurations. Furthermore, iTollerud et al.l ()2011f ) find 
that, compared to other R-band selected LMC-MW 
analogs in SDSS, the LMC is unusually blue in color. 
This fact is reconcilable in a first infall scenario where 
the LMC is accreted recently, because it has been able 
to retain enough gas to maintain a high star forma- 
tion rate today (see also iWetzel et al.l 120121 ) . However , 
ITollerud et 311(12011!) also find that the LMC's color is 
unusual c ompared to R-ban d selected LMC analogs in 
the field. iBesla et al.l ()2012D present a theory in which 
the star formation rate of both the LMC and SMC has 
increased recently owing to interactions between these 
two galaxies. Given that LMC-SMC configurations are 
rare in general, it is thus unsurprising that the overall 
star formation rate of the LMC is higher than average 
isolated analogs. 

Even in the case for the most long-lived LMC-MW con- 
figuration presented here (i.e., low LMC mass and high 
MW mass) the LMC makes at best ^ 2 past pericentric 
passages (see Figure [H]). The long implied orbital periods 
excee d the lifetime of the Magellanic Stream (1-2 Gyr ) 
(e.g.. lHeitsch fc Putmanll2009l: iKeres fc Hernauistll2009D . 
This conclusion is consistent with our arguments for the 
need for a new formation mechanism for the Magellanic 
Stream that does not rely on a previous pericentric ap- 
proach about the MW. Moreover, such a combination of 
LMC and MW masses make it difficult for the LMC and 
SMC to have been a long-lived binary system. 

The existence of the LMC-SMC system as a binary for 
more than 2 Gyr in the past favors a combination of a 
high mass LMC and an intermediate or low mass MW. 
Our work has found that such configurations always yield 
first infall orbital configurations. We show an example 
of an orbit that fulfills all the criteria set out above in 
Figure 1131 This shows that viable orbital solutions can 
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be found that are consistent with both the HST PMs 
and cosmological expectations. The new SMC PM now 
ahows the past orbit of the SMC to cross-over the lo- 
cation of the HI in the Stream, which was very difficult 
to obtain with the old K2 PMs and is consistent with a 
picture in which the Stream is formed pr imarily by the 
removal of material from the SMC (GN96:lConnors et alj 
[20061: IBesla et alJl20ia[20TllDiaz fc Bekkill201lL I20m 

7.3. Future Work 

Recently the sample of QSOs behind the Magellanic 
Clouds has increased drastically. There are now well 
over 200 newly identified QSOs behind the LMC and 29 
behind the SMC. The new QS Os were selected primar- 
ily ba sed on their mid-IR colors (Kozlowski fc KochaneU 
l2009f) . The majority of them were spectroscopically con- 
firmed from follow-up data tak en with 2DF/AA0MEGA 
(jKozlowski et al.l 120 111 l2012t ). while a smaller subset 
(around 40) were confirmed with Magellan/IMACS (by 
N. Kallivayalil fc M. Geha). PM studies of fields around 
a carefully chosen subset of these QSOs could vastly im- 
prove our knowledge of the structure and dynamics of 
the Clouds. 



With a larger sample of QSOs homogeneously dis- 
tributed behind the SMC, we can investigate its rotation 
and structure, something we have been unable to tackle 
with only 5 QSO fields thus far (and with only 3 QSOs 
currently imaged with WFC3/UVIS). A larger number of 
QSOs will also improve our knowledge of its COM PM. 

For the LMC, a large number of fields with PMs would 
allow us to better sample the rotation curve, including 
its inner slope. The increased coverage would also al- 
low u s to combine our PMs with radial velocity studies 
(e.g.. iParisi et al.ll2010l : lOlsen et al. 2011), thereby con- 
stituting a unique 3D data-set. With a large number of 
QSO fields we may also be able to distinguish between 
the kinematics of different stellar populations within the 
LMC directly. 

Finally, with exposures of sufficient depth, it is possi- 
ble to m easure abso lute PMs with resolved background 
galaxies (jSohn et al., 2012). so we should not be confined 
to fields with QSOs. We are in the process of identifying 
fields in the Clouds with deep HST/ ACS exposures in an 
earlier epoch, which should have numerous background 
galaxies that can be used to yield an accurate absolute 
PM of the stars in the field. 



APPENDIX 
WFPC2/PC DATA 

A third epoch of HST imaging of our QSO fields was originally due to execute on ACS/HRC in the period July 2007 
to November 2008. However, due to the failure of ACS at that time, the observations were executed on the Planetary 
Camera (PC) of the WFPC2 instrument (PID 11201). All of the 40 targets (34 LMC & 6 SMC QSO fields) that had 
been observed with ACS/HRC in epoch 1 were observed with WFPC2 in snapshot mode, providing a typical time 
baseline of 5 years. The observational approach was similar to that in epochs 1 fc 2. We used the l/-band (F606W) 
filter and a 5 or 6-point dither pattern. Total exposure times ranged from 2.8 to 20 minutes. 

The data were analyzed with similar techniques as used for the first two ACS epochs. A discussion of the analysis 
and results for a subset of the fields was presented in iKallivavalil et all ()2009b[) . In Figure US] we show the position 
of the quasar over time, in a reference frame based on the 1st epoch images, for 4 randomly chosen LMC fields. The 
scatter among positions within the same epoch shows that the RMS error in the position of the quasar is roughly 3 
times as large for WFPC2 as for ACS (~ 0.021 HRC pixels versus 0.008 HRC pixels, respectively). The lower accuracy 
for WFPC2 is due to a combination of several factors. These include the larger pixel size, lower sensitivity, and less 
well calibrated and stable geometric distortions for WFPC2 compared to ACS/HRC. 

The 40 fields for which WFPC2 data were obtained fall in three categories: 

(1) 14 Fields for which two epochs of ACS data exist, but no WFC3 data. For these fields the WFPC2 data provide 
a useful consistency check, as evident from Figure [TH The WFPC2 data more than double the time baseline and fall 
on the extrapolation of the ACS measured PMs. The use of a different i nstrument and the additi on of a third epoch 
rule out a variety of potential systematic errors. However, as discussed in IKallivavalil et "aD (|2009b( ). the WFPC2 data 
help very little to reduce the random PM errors in each field. 

(2) 13 Fields for which ACS data exist, as well as one epoch of WFC3 data. For these fields the WFC3 data provide 
superior accuracy and a longer time baseline. Hence, the WFPC2 data do not provide anything in terms of additional 
insights or reduced PM errors. 

(3) 13 Fields for which only a first epoch of ACS data exists. For 11 of these fields the data provide the opportunity 
to determine a PM where none yet exists (2 fields are unsuitable as discussed in Appendix B). This might in principle 
help to obtain better spatial coverage of measurements within the LMC and SMC. However, we have found that the 
PM accuracy that can be achieved for these fields is relatively low. The random PM errors for 5-year ACS-WFPC2 
measurements are a factor 2-4 higher than for 2-year ACS-ACS measurements. 

These considerations indicate that the WFPC2 data are of insufficient quality to help better constrain the analysis 
presented in this paper. For this reason, we have not attempted to analyze the WFPC2 data to the same level of detail 
and completion as for the other instruments, and we do not report and use the WFPC2 results in this paper. 

OBSERVED FIELDS UNSUITABLE FOR PM DETERMINATION 

We found that two QSO fields observed with WFC3 proved unsuitable for PM determination of the Magellanic 
Clouds. These fields were not part of the Kl or K2 sample, but only had one epoch of previous ACS data (which 
had never been fully analyzed). We briefly discuss here the reasons for excluding these fields from the sample for the 
present paper. 
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Figure 15. Positions, in a reference frame based on the 1st epoch images, of the quasar over time for 4 randomly chosen LMC fields 
showing ACS epoch 1 (crosses), ACS epoch 2 (triangles) and WFPC2/PC (open squares). The scatter per epoch gives the relative size of 
the RMS errors in the position of the quasar relative to the star-field, given that the observation has been repeated many times per epoch. 
Each panel is 1 X 1 pixel (1 pix = 28 mas). North is in a different direction in each panel. 

RX J0524.0-7011: This ROSAT X-Ray source at (a, 5) = (05fe2 4m02.3s, ~70dllm09s) was optically identified 
as an AGN at z = 0.15 in th e LMC background by ICowlev et al.l (17984). It was cross- identified as MACHO ID 
006.07059.0207 by Ge ha et all ()2003D . It has the lowest redshift in our sample, and is also relatively bright, with 
V = 17.7. In our WFC3 images the AGN is clearly resolved, and a spiral host galaxy is visible surrounding a point-like 
nucleus. The astrometric techniques used in this paper are not designed for extended sources, and application of these 
techniques was found not to produce reliable PM results for this target. We therefore excluded it from the sample. 
RX JOlll. 7-7250: This ROSAT X-Ray source was included by Tinney et al. (1997) in a study aimed at finding 
QSOs behind the SMC. They aligned an optical spectrograph slit along the two bluest objects within a ^ 20" circle 
(believed to be a good guess for the ROSAT positional error) around the source. Neither of the objects turned out 
to be a background AGN. However, they detected a serendipitous background emission-line object along the slit at 
redshift z = 0.197. The object was called QJOlll-7249, and identified as an optical QSO. This seems tentative at 
best, given that the spectrum showed only Balmer emission lines (which leaves open the possibility for alternative 
categorizations). The emission-line object was 23" from the ROSAT position at the time, only marginally consistent 
within the uncertainties. Subsequently, the X-ra y source was inc luded as #348 in the ROSAT PSPC catalogue of 
X-ray sources in the SMC region compiled by H aberl et al.l ()200CI( ). They provided an improved position of {a,S) = 
(01/illm40.9s — 72c?50to28s) (J2000.0), with a 90% confidence random uncertainty of 4.1" and an additional systematic 
uncertainty of ^ 7". This position is 41" from QJOlll-7249, and the sources therefore are not associated. Nonetheless, 
most literature has continued to treat the sources as one and the same. Our HST observations were centered on 
QJOlll-7249. Hence, the actual X-ray source did not fall in the ACS/HRC field of view. Our PM analysis of the 
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field did not reveal any point-like source with a reflex motion typical of a background reference source. Therefore, we 
conclude that QJOl 11-7249 is not suitable for SM C PM determination, and exclude it from the sample. Interestingly, 
this same source was used bv lCosta et al.l ()2011[ ) in their ground-based study of the SMC PM. They obtained a very 
discrepant motion for this SMC field, {hw,IJ'n) = (—0.06, —0.71) mas yr~^ , com pared to the i r aver age based on the 
other fields, (—0.93, —1.25) mas yr~^ . This supports our interpretation, although iCosta et al.l (|2011| ) argue it may be 
due to streaming motions in the SMC. 

Both of the targets discussed in this appendix are atypical in our sample, given that they were selected a s previously 
known X-ray sources. The large majority of our sample studied in this paper is composed of QSOs from iGeha et al.l 
(|2003D . which were identified from their variability characteristics in the MACHO database, with subsequent spectro- 
scopic follow-up. 
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Telescope Science Institute (STScI), which is operated by the Association of Universities for Research in Astronomy, 
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Table 1 

LMC Observations and PM Results 



ID RA DEC T Mime N PM of FIELD AS OBSERVED LMC PM(CM) estimate 





(H,M,S) 


(deg, ' 




(min) 


(vr) 






2-epoch PMs 

{ mas yr^-"- ) 






3-epoch PMs 
( mas yr~^ ) 




( 


mas yr^"*" ) 


LI 


5 


47 50.2 


-67 


28 1.3 


2.6 


7.1 


8 


-1 


.491 


0.689 


0.076 





.069 


-1.604 


0.714 


0.029 





.032 


-1 


.820 


0.209 


L2 


4 


46 11.1 


-72 


5 9.0 


6.6 


7.7 


13 


-2 


.163 


-0.246 


0.061 





.054 


-2.051 


-0.235 


0.018 





.027 


-1 


.910 


0.181 


L3 


5 


16 28.9 


-68 


37 


1.8 


6.6 


7.1 


45 


-1 


.827 


0.391 


0.047 





.042 


-1.798 


0.320 


0.014 





.016 


-1 


.924 


0.252 


L4 


5 


16 36.8 


-66 


34 


35.8 


12.6 


7.6 


9 


-1 


.781 


0.196 


0.083 





.093 


-1.734 


0.147 


0.048 





.030 


-1 


.974 


0.100 


L5 


5 


15 36.1 


-70 


54 


0.8 


- 


1.9 


41 


-2 


.008 


0.112 


0.079 





.079 


- 










-1 


.767 


0.143 


L6 


6 


2 34.3 


-68 


30 


41.1 


13.3 


7.2 


8 


-1 


.750 


0.897 


0.158 





.177 


-1.664 


1.058 


0.045 





.055 


-1 


.873 


0.392 


L7 


5 


16 26.3 


-69 


48 


19.0 


- 


1.8 


127 


-1, 


.942 


0.226 


0.053 


0, 


.053 


- 










-1, 


.824 


0.207 


L8 


5 


54.0 


-66 


43 


59.8 


- 


2.0 


39 


-1, 


.883 


-0.0232 


0.062 


0, 


.062 


- 










-2, 


.068 


0.247 


L9 


4 


53 56.5 


-69 


40 


35.4 


- 


2.0 


73 


-1, 


.985 


-0.058 


0.079 


0, 


.079 


- 










-1, 


.892 


0.394 


LIO 


5 


49 41.6 


-69 


44 


15.1 




1.9 


65 


-1, 


.704 


0.648 


0.087 


0, 


.087 












-1, 


.802 


0.062 


Lll 


5 


57 22.4 


-67 


13 


21.5 


2.9 


6.6 


8 


-1, 


.591 


0.895 


0.174 


0, 


.078 


-1.664 


0.814 


0.041 


0, 


.031 


-1, 


.899 


0.233 


L12 


5 


1 46.7 


-67 


32 


39.8 


13.3 


6.4 


16 


-1, 


.903 


-0.201 


0.083 


0, 


.121 


-1.817 


-0.065 


0.033 


0, 


.020 


-1, 


.980 


0.231 


L13 


5 


20 57.0 


-70 


24 


52.6 




1.5 


85 


-2, 


.199 


0.543 


0.093 


0, 


.093 












-1, 


.980 


0.409 


L14 


5 


2 53.7 


-67 


25 


45.0 


13.3 


7.1 


18 


-1, 


.664 


-0.091 


0.108 


0, 


.140 


-1.707 


-0.054 


0.034 


0, 


.031 


-1, 


.885 


0.212 


L15 


5 


14 12.1 


-70 


20 


25.8 




1.1 


75 


-2, 


.512 


0.553 


0.127 


0, 


.127 












-2, 


.273 


0.628 


L16 


5 


4 36.0 


-66 


24 


15.7 


10.7 


7.1 


13 


-1, 


.561 


0.166 


0.063 


0, 


.077 


-1.684 


0.079 


0.031 


0, 


.035 


-1, 


.890 


0.270 


L17 


5 


25 14.4 


-65 


54 


45.7 




1.2 


17 


-1, 


.773 


1.247 


0.290 


0, 


.290 












-2, 


.009 


1.081 


LIS 


5 


30 26.8 


-66 


48 


52.9 




1.6 


22 


-1, 


.767 


0.716 


0.150 


0, 


.150 












-2, 


.001 


0.434 


L19 


5 


31 59.7 


-69 


19 


51.5 




1.1 


114 


-1, 


.636 


0.780 


0.150 


0, 


.150 












-1, 


.706 


0.340 


L20 


4 


56 14.3 


-67 


39 


9.0 




2.5 


48 


-1, 


.826 


-0.139 


0.057 


0, 


.057 












-1, 


.939 


0.260 


L21 


5 


10 32.5 


-69 


27 


15.5 




1.1 


129 


-2, 


.042 


-0.113 


0.118 


0, 


.118 












-1, 


.964 


0.043 


L22 


5 


20 56.5 


-65 


39 


4.8 


17.7 


6.9 


23 














-1.681 


0.298 


0.020 


0, 


.036 


-1, 


.914 


0.206 



Note. — The LMC field number, using the same notation as in Kl, and RA/DEC for the QSOs (eolumns 1, 2 & 3). Column 4 lists the total exposure times 
in minutes for e3. Column 5 lists the time baseline, in years, between el and e3. Column 6 lists the number of stars used in the final transformations after all 
cuts and iterations have been applied. Columns 7—10 list the observed PMs and errors for the ACS reanalysis (see ^2.3| l . Columns 11—14 list the corresponding 
observed PMs and errors obtained from the 3-epoch analysis (see ^2.5| l . Columns 15 &z 16 list the final COM PM estimates from each QSO field obtained from 
the procedure described in i[3?2l Note that L22 is a new field, with only el and e3 data. 



Table 2 

SMC Observations and PM Results 



ID RA DEC T Atime N PM of FIELD AS OBSERVED SMC PM(CM) estimate 



















(5/ijv 






(5/jjv 




















2-epoch PMs 






3-epocli PMs 










(H,M,S) 


(deg, ") 


(min) 


(yr) 






( mas yr~i ) 






( mas yr"'- ) 




( 


mas yr"-"- ) 


SI 


00 51 17.0 


-72 16 51.3 




1.9 


42 


-0.682 


-1.288 0.100 


0.100 








-0.738 


-1.244 


S2 


00 55 34.7 


-72 28 33.9 


17.7 


7.6 


25 


-0.748 


-1.246 0.070 


0.068 


-0.722 


-1.214 0.032 


0.024 


-0.760 


-1.185 


S3 


01 02 14.5 


-73 16 26.6 


17.7 


7.7 


36 


-0.893 


-1.397 0.097 


0.101 


-0.679 


-0.974 0.026 


0.028 


-0.691 


-0.973 


S4 


00 36 39.7 


-72 27 42.0 




2.8 


10 


-0.460 


-1.114 0.109 


0.109 








-0.579 


-1.031 


S5 


01 02 34.7 


-72 54 23.8 


13.3 


6.8 


30 


-1.046 


-1.072 0.084 


0.083 


-0.806 


-1.199 0.017 


0.038 


-0.815 


-1.195 



Note. — 



Same columns as in Table fl] but for the SMC. 
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Table 3 

Kinomatical determinations of LMC geometry and rotation parameters 







Faper 11 


vuMUz 


Olsenll 


(1) 


(2) 


(3) 


(4) 


(5) 


inc 


deg 


39.6 ±4.5 


34.7 ±6.211' 


34.7 ± 6. 2W 


theta 


deg 


147.4 ± 10.0 


129.9 ± 6.0 


142 ±5 


RA 


deg 


78.76 ±0.52 


81.91 ± 0.98 


81.91 ± 0.98121 


DEC 


deg 


-69.19 ±0.25 


-69.87 ±0.41 


-69.87 ±0.4lPl 




kpc 


1.18 ±0.48 


4.0 ± 0.3(51 


2.4 ±0.1 




km s^-"^ 


76.1 ±7.6M 


49 ± 2 [5] 


87 ± 5M 


^sys 


km s^-*^ 


262.2 ± 3. 4^1 


262.2 ± 3.4 


263 ±2 


m - M 


mag 


18.50 ±0.10M 


18.50 ±0.10M 


18.50 ±0.10M 



Note. — Column (1) lists the following quantities: inclination of the LMC 
disk plane; position angle of the line of nodes; (RA,DEC) of the rotation center; 
turnover radius Ro and asymptotic velocity amplitude Vo of the rotation curve; 
systemic LOS velocity, Vsys', and distance modulus, m — M. Column (2) lists 
the corresponding units. Column (3) lists the values inferred from the new PM 
data as described in Paper IL Columns (4) and (5) li s t the values inferred from 
the LOS velocity studies of vdM02 and ioisen et all I|2011D . 

[^1 Value from Ivan de r Marcl fc Cioni' f2001]) , used but not independently de- 
termined by vdM02 or Olscn ot al. (2011.). 

PI Val ue from vdM02, used but not independently determined bv lOlsen et al.l 

1^1 The rotation curve is parameterized so that it rises linearly to velocity Vo 
at radius _Ro, and then stays flat at larger radii. 

Applies to a mix of stellar populations (each with different asymmetric drift) 
as shown in Figure 6 of Kl. 

1^1 Applies to (old) carbon stars (with large asymmetric drift). Determined 
from Table 2 of vdM02. Used an LMC COM PM value that pre-dates the now 
current HST values. 

[^1 Applies to (young) red supergiants (with little asymmetric drift). Used the 
LMC COM PM from P08. 

PI Value from vdM02, used but not i ndependently determined in Paper II. 

Value from IFreedman et al.l l|2001f ). used but not independently determined 
from kinematics. 



Table 4 

Summary of Recent LMC and SMC PM Measurements 





LMC 


SMC 
















Work 


( mas yr~i ) 


( mas yr"-*- ) 


( mas yr~^ ) 


( mas yr"-*- ) 


Data 


This Paper 


-1.910 ± 0.020 


0.229 ± 0.047 


-0.772 ± 0.063 


-1.117 ± 0.061 


HST 3-epoch 


(vdM02 model) 


-1.899 ±0.017 


0.416 ±0.017 






HST 3-epoch 


Kl & K2 


-2.03 ±0.08 


0.44 ± 0.05 


-1.16 ±0.18 


-1.17 ±0.18 


HST 2-epoch 


P08 


-1.956 ±0.036 


0.435 ± 0.036 


-0.754 ±0.061 


-1.252 ± 0.058 


HST 2-epoch 


Costa et al. 


-1.72 ±0.13 


0.50 ±0.15 


-0.93 ±0.14 


-1.25 ± 0.11 


2.5m du Pont 


Vieira et al. (2010) 


-1.89 ±0.27 


0.39 ±0.27 


-0.98 ±0.30 


-1.10 ±0.29 


SPM 



Note. — PMs of the LMC (columns 3 and 4) and SMC (columns 5 and 6) COM. Column (1) indicates the source 
of the result, and column (6) the type of data that was used. The first line (labeled "This Paper") is the final result 
from the present paper, which uses the fit to the PM rotation field from Paper II. Ours is the first study that propagates 
the uncertainties in the geometry, center and rotation of each Cloud into the COM PM estimates. For this reason, other 
studies have generally underestimated their error bars. The second line shows the LMC COM PM estimate that is obtained 
when the LMC orientation and center are are kept fixed to the same vdM02 values that were used by most other authors 
(see Section |4.1| l. For Costa et al. we list for the LMC the average of their 2009 results for V^ot — 50 km/s and 120 km/s, 
respectively; for the SMC we list their 2011 result. 
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Table 5 

Galactocentric Velocities from HST measurements 



Line 


Galaxy 


solar 


PM 


vx 


vy 


vz 


Vlot 


frad 


ft an 










km/s 


km/s 


km/s 


km/s 


km/s 


km/s 


(1) 


LMC 


new 


3-epoch 


-57 ± 13 


-226 it 15 


221 it 19 


321 it 24 


64 it 7 


314 it 24 


\ ) 


LMC 


lAU 


3-epocli 


c;q -I- 1 9 


9f^9 -1- 1 ^ 


221 it 19 


'\AC\ A- 9*? 
o^u in 


ou in o 




\ > 


LMC 




pnnrh vrlMn9 

O cUO^jil Vt_liVlVJ^ 


— 77 it 8 


—224 it 14 


227 it 18 


OiO ^ zo 


65 it 5 


Q99 -L 94 
ozz m zt 




LMC 


lAU 




— 78 it 8 


—250 it 14 


227 it 18 


347 it 23 


87 it 5 




(5) 


LMC 


lAU 


P08 


-83 it 11 


-258 it 15 


238 it 20 


361 it 25 


88 it 5 


350 it 26 


u 


LMC 


lAU 


Kl 


-86 it 14 


-268 it 18 


252 it 25 


378 it 31 


89 it 5 


367 it 31 


(7) 


SMC 


new 


3-epoch 


19 it 18 


-153 it 21 


153 it 17 


217 it 26 


-11 it 5 


217 it 26 


(8) 


SMC 


lAU 


3-epoch 


18 it 17 


-179 it 21 


153 it 17 


236 it 26 


6it4 


236 it 26 


9) 


SMC 


lAU 


P08 


23 it 16 


-197 it 22 


166 it 17 


259 it 26 


7it4 


259 it 26 


(10) 


SMC 


lAU 


K2 


-86 it 49 


-248 it 46 


150 it 39 


302 it 57 


23 it 7 


301 it 57 


(11) 


SMC-LMC 




3-epoch 


76 it 22 


73 it 26 


-68 it 25 


128 it 32 


112 it 32 


61 it 16 


(12) 


SMC-LMC 




3-epoch vdM02 


81 it 19 


73 it 26 


-89 it 25 


143 it 31 


134 it 32 


50 it 15 


(13) 


SMC-LMC 




P08 


106 it 20 


61 it 26 


-72 it 26 


145 it 30 


136 it 33 


46 it 18 


(14) 


SMC-LMC 




K2 


OitSl 


20 it 49 


-103 it 46 


127 it 46 


77 it 50 


88 it 45 



Note. — The lines Cl)-(6) list the LMC velocity, hnes (7)-{10) the SMC velocity, and hnes (11)-(14) the relative velocity of the 
SMC with respect to the LMC. Column (1) lists a line identifier. Column (2) lists the galaxy name. Column (3) lists the source of 
the velocities used to correct for solar refl^ex motion: "lAU" uses the lAU value of Vb ^ 220 km s"^ and the lDehnen & Binncv ( 1998) 
solar peculiar velocity; "new" uses the improved iMcMillanI ||201in value of Vo — 239 ± 5 km and the improved [Schonrich ct_alj 

|||2 01C|) solar peculiar velocity. Column (4) lists the assumed PM value, taken from the list of observations in Table [4] "3-epoch" 
corresponds to the line labeled "This Paper" (which uses the new data with the LMC geometry fit from Paper II, with uncertainties 
included); "3-epoch vdM02" corresponds to the line labeled "(vdM02 model)" (which uses the new data with the fixed geometry 
parameters from vdM02). Columns (5)-(7) list the Galactocentric velocity coordinates {vx , vy , vz) ■ Columns (8)-(10) lists the 
total length of the velocity vector, the radial component, and the transverse component, respectively. Uncertainties were calculated 
using a Montc-Carlo scheme that propagates all relevant uncertainties in the position and velocity of both the Clouds and the Sun. 
Distance uncertainties are based on Am — M — 0.1. Velocity uncertainties in the Galactocentric frame are highly correlated, because 
uncertainties in the LOS direction are much smaller than in the transverse direction. 



Table 6 

MW and LMC Properties 



Galaxy 


Total Mass 

[Mq] 




-Rvir 

[kpc] 


A^disk 


Plummer Softening 
[kpc] 


MW 


1 X 10^2 


9.86 


261 


6.5 X 10"^ 




MW 


1.5 X 10^2 


9.56 


299 


5.5 X 10^° 




MW 


2 X 10^2 


9.36 


329 


5.0 X lOio 




LMC 


3 X 10^° 








8 


LMC 


5 X 10^'' 








11 


LMC 


8 X 10^'' 








14 


LMC 


1 X 10" 








15 


LMC 


1.8 X 10" 








20 


LMC 


2.5 X 10" 








22.5 



Note. — Different MW and LMC mass models used in the orbital calcula- 
tions. Cvir is the halo concentration, and i?vir is the virial radius of the MW 
in each case. The mass of the MW bulge is kept fixed at 1 X 10^° M© and 
the Hernquist scale radius for the MW bulge is kept fixed at 0.7 kpc. The 
MW's exponential disk scale radius is kept fixed at 3.5 kpc. The mass of the 
MW disk is varied to get the observed circular velocity at the solar circle as 
discussed in the text. 



